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Abstract 

The  nitrogen- vacancy  (NV)  center  in  diamond  has  emerged  as  a  promising  platform  for  high- 
sensitivity,  vector  magnetic  field  detection  and  high  spatial  resolution  magnetic-field  imaging 
due  to  its  unique  combination  of  optical  and  spin  properties.  NV  diamond  magnetometry  has 
enabled  a  wide  array  of  applications  from  the  noninvasive  measurement  of  a  single  neuron 
action  potential  to  the  mapping  //T-fields  in  ji m-size  meteorite  grains.  To  further  improve  the 
magnetic  sensitivity  of  an  ensemble  NV  magnetometer,  the  growth  and  processing  of  the  host 
diamond  must  be  taken  into  account.  This  thesis  presents  a  systematic  study  of  the  effects  of 
diamond  processing  on  bulk  chemical-vapor-deposition  diamond.  In  particular,  NV  charge- 
state  composition  and  spin  decoherence  times  are  measured  for  diamonds  irradiated  with  1 
MeV  electrons  at  doses  of  1  x  1015  —  5  x  1019  e“/cm2  and  thermally  annealed  at  temperatures 
of  850°C  and  1250°C.  The  study  provides  an  optimal  range  for  diamond  processing  and  shows 
the  quenching  of  the  NV  center  at  high  irradiation  dosage  from  the  creation  of  additional 
vacancy-related  defects. 

Thesis  Supervisor:  Paola  Cappellaro 

Title:  Esther  and  Harold  E.  Edgerton  Associate  Professor  of  Nuclear  Science  and  Engineering 
Thesis  Supervisor:  Danielle  A.  Braje 

Title:  Assistant  Group  Leader,  Quantum  Information  and  Integrated  Nanosystems 
MIT  Lincoln  Laboratory 


3 


4 


Acknowledgments 


My  time  at  MIT  has  been  two  of  the  most  challenging  years  of  my  life,  pushing  my  mental 
stamina  far  beyond  anything  I  could  have  imagined.  The  support  of  many  individuals,  both 
in  the  local  area  and  across  the  country,  are  the  main  reason  for  my  accomplishments  while 
pursuing  this  degree. 

All  of  this  would  not  have  been  possible  in  the  first  place  were  it  not  for  the  support  of 
my  sponsors.  I  am  extremely  thankful  for  the  United  States  Air  Force  for  allowing  me  to 
pursue  graduate-level  education  for  my  first  assignment  as  a  commissioned  officer.  I  hope  to 
use  the  many  skills  I  have  learned  here  to  assist  my  future  colleagues  when  I  begin  my  work 
as  a  physicist  for  the  Air  Force  Research  Laboratory  in  New  Mexico.  I  am  also  thankful  for 
the  National  Science  Foundation  and  the  MIT  Lincoln  Laboratory  for  funding  my  graduate 
education.  The  material  in  this  thesis  is  based  upon  work  supported  by  the  National  Science 
Foundation  Graduate  Research  Fellowship  under  Grant  No.  1122374. 

Words  alone  cannot  express  how  thankful  I  am  for  my  two  advisors,  Dr.  Danielle  Braje  of 
Lincoln  Laboratory  and  Professor  Paola  Cappellaro  of  MIT.  Having  near-zero  experimental 
experience  when  I  started,  they  patiently  guided  me  through  how  to  take  and  present  data, 
in  addition  to  valuable  life  lessons  on  what  being  an  experimental  physicists  is  like.  I  look 
forward  to  continued  mentorship  from  them. 

I  am  highly  grateful  for  the  two  technical  staff  members  of  Lincoln  I  frequently  worked 
with,  Dr.  Linh  My  Pham  and  Dr.  John  Barry.  Linh  answered  the  many  optics  questions 
I  had  in  my  first  year  and  taught  me  everything  I  know  about  characterizing  NV  centers 
in  my  second  year.  If  I  ever  had  a  problem  with  the  experiment  that  despite  hours  of 
troubleshooting  I  could  not  fix,  Linh  always  knew  the  right  fix  and  would  explain  it  in 
a  clear  and  pedagogical  way.  I  will  also  miss  hearing  her  insightful  critiques  on  several 
movies,  especially  those  from  the  DC  and  Marvel  universe.  I  had  the  pleasure  of  working 
with  John  my  second  year,  especially  towards  the  end  as  we  attempted  to  figure  out  EPR 
measurements  and  spend  several  hours  at  MIT  on  the  weekends.  Watching  him  troubleshoot 
the  EPR  machine  and  spectrometer  taught  me  valuable  lessons  on  how  to  take  data  efficiently 
and  make  sure  all  possible  instrumentation  errors  are  accounted  for.  Most  importantly  I  am 


5 


thankful  for  their  constant  advice  on  how  to  be  a  professional  experimental  physicist,  from 
teaching  me  to  focus  my  attention  on  what  is  truly  important  for  my  experiment  to  always 
challenging  me  to  consider  all  ways  data  acquisition/processing  can  be  made  more  efficient. 

I  am  thankful  for  sharing  a  lab  bench  with  Erik  Eisenach  and  one  bench  over  with 
Christopher  McNally.  The  constant  stream  of  laughter  and  group  lunches  made  every  day 
at  Lincoln  worth  enjoying.  I  wish  Erik  the  best  as  he  continues  to  get  a  doctorate  at  MIT 
and  for  Chris’s  future  graduate  school  endeavors.  The  “graduates”  of  the  Quantum  Sensing 
Laboratory  I  am  also  thankful  for;  the  assistance  of  Hannah  Clevenson,  Matt  Trusheim,  and 
Edward  Chen  were  invaluable  during  my  first  year  at  Lincoln.  I’m  especially  fond  of  Ed’s 
intensity  towards  physics,  epitomized  by  his  quote  “you  are  always  at  war  with  science”. 

Shifting  to  the  MIT  campus,  I  am  unquestionably  grateful  for  the  time  spent  in  the 
Quantum  Engineering  Group  (QEG).  I  have  fond  memories  of  the  late-night  problem  set 
parties  of  22.51  and  8.421  with  Calvin  Sun  and  Yixiang  Liu.  I’m  also  thankful  for  the 
assistance  of  David  Layden  in  8.511  and  8.422  during  my  second  year  of  classes.  The  insight 
and  welcoming  atmosphere  provided  by  Mo  Chen,  Akira  Sone,  Ken  Xuan  Wei,  Kasturi  Saha, 
“professor”  Alex  Cooper-Roy,  Ashok  Ajoy,  and  Luca  Marseglia  made  every  interaction  with 
QEG  something  to  look  forward  too.  I’m  excited  to  see  what  research  they  will  accomplish 
in  the  future. 

Outside  of  QEG  are  several  other  fellow  students  I  enjoyed  spending  time  with.  There’s 
no  one  else  I’d  rather  go  through  the  course  22  modules  with  than  Chris  Wink  and  Malik 
Wagih.  Those  also  involved  lots  of  late  nights  trying  to  understand  just  what  was  really 
being  asked  in  problem  sets.  I’m  also  happy  to  have  met  Jules  Stuart  (“153”)  of  the  ion¬ 
trapping  group;  I  will  miss  his  humor  from  working  on  solid-state  class  work  and  on  MIT 
in  general.  Finally,  I  had  the  tremendous  pleasure  of  taking  both  of  Professor  Wolfgang 
Ketterle’s  atomic  and  optical  physics  courses;  from  these  I  have  gained  an  appreciation 
for  AMO  science  and  a  baseline  knowledge  for  my  next  assignment  when  I  will  work  with 
AMO  experiments.  Professor  Isaac  Chuang,  who  co-managed  the  courses,  introduced  me  to 
the  flipped  classroom  paradigm  with  half  of  the  course  being  on  edX,  a  new  and  insightful 
learning  experience  I  will  not  forget. 

I’ve  also  been  privileged  to  interact  with  the  Air  Force  scientific  community.  I  am  end- 


6 


lessly  grateful  for  the  continued  mentorship  of  my  undergraduate  advisor  Lt  Col  Mario  Serna, 
who  helped  me  keep  the  bigger  picture  in  mind  with  regards  to  MIT  and  the  Air  Force.  I 
had  the  privilege  of  sharing  my  research  with  a  few  Air  Force  Research  Laboratory  sites; 
I  am  thankful  for  being  allowed  to  present  my  research  to  Dr.  Luke  Bissell  and  Dr.  Paul 
Alsing.  I’m  also  honored  that  Dr.  Spencer  Olson  and  Lt  Col  Serna  were  able  to  listen  in  to 
my  final  thesis  presentation  from  abroad. 

Lastly,  I  am  thankful  for  the  wonderful  friendships  I  have  made  at  the  MIT  Baptist 
Student  Fellowship  and  the  continued  friendships  from  my  time  at  the  United  States  Air 
Force  Academy.  And  of  course,  I  am  extremely  thankful  to  my  family  back  home  in  Florida, 
to  now  include  my  brother’s  fiance.  They  always  encourage  me  to  try  harder  and  redefine 
what  I  consider  my  limits  to  be.  In  return,  I  hope  I  was  able  to  convince  them  of  the 
usefulness  of  the  nitrogen-vacancy  center! 


7 


Disclaimers 


The  research  presented  in  this  thesis  is  jointly  supported  by  the  National  Science  Foundation 
and  MIT  Lincoln  Laboratory. 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or  the  U.S.  Gov¬ 
ernment. 

This  material  is  based  upon  work  supported  by  the  Assistant  Secretary  of  Defense  for 
Research  and  Engineering  under  Air  Force  Contract  No.  FA8721-05-C-0002  and/or  FA8702- 
15-D-0001.  Any  opinions,  findings,  conclusions  or  recommendations  expressed  in  this  mate¬ 
rial  are  those  of  the  author  and  do  not  necessarily  reflect  the  views  of  the  Assistant  Secretary 
of  Defense  for  Research  and  Engineering. 

Any  opinion,  findings,  and  conclusions  or  recommendations  expressed  in  this  material 
are  those  of  the  author  and  do  not  necessarily  reflect  the  views  of  the  National  Science  Foun¬ 
dation. 


This  thesis  is  dedicated  to  my  parents,  Mark  and  Tracey 
Alsid,  and  to  my  brother  Brett.  I  am  forever  thankful 
for  their  love  and  support. 


9 


10 


Contents 


1  Magnetometry  and  the  Nitrogen- Vacancy  Center 

11.1  Introduction! . 

1.2  Nitrogen- Vacancy  Center  Physicsj . 

11.2.1  Electronic  Structure! . 

1.2.2  Photoluminescence  Spectra] . 

1.2.3  The  Neutrally-Charged  NY  Center . 

1.2.4  Spin-State  Dependent  Fluorescence  and  ODMR 

1.3  NY  Magnetometry| . 

1.3.1  Magnetic  Sensitivity . 

1.3.2  Optimizing  Sensitivity . 


2  Diamond  Processing! 

12.1  Overviewl . 

2.2  Diamond  Processing! . 

2.2.1  Synthetic  Diamond  Generation  .  .  . 

2.2.2  Nitrogen  Incorporation . 

2.2.3  Vacancy  Production  . 

2.2.4  Thermal  Annealing) . 

12.3  Diamonds  used  in  this  Thesis! . 

12.4  Conclusion! . 

3  Charge  State  Fluorescence  Characterization 

13.1  Overviewl . 


25 

25 

26 

28 

29 

31 

33 

36 

37 

37 

41 

41 

43 

43 

45 

46 

47 

47 

50 

51 

51 


11 


52 


3.2  Charge  State  Dynamics  and  Magnetometry . 

13.3  Measurement,  Procedurel . 

3.4  Charge  State  Fluorescence  Measurements . 

13.4.1  Total  Fluorescence  vs.  Irradiation  Dosel . 

13.4.2  Fluorescence  Modell . 

3.4.3  Total  Fluorescence  vs.  Optical  Power . 

3.4.4  Charge  State  Fluorescence  Ratio  vs.  Irradiation  Dose 

3.4.5  Charge  State  Fluorescence  Ratio  vs.  Optical  Power  . 

13.5  Con.dnsi.onl . 

14  Coherence  Time  Characterization! 

14.1  Qverviewl . 

4.2  Coherence  Time  and  Magnetometry] . 

4.3  Spin-Dephasing  Time  T2*  Measurements . 

14.3.1  Measurement  Procedurel . 

14.3.2  Results  and  Limitations! . 

4.4  Spin-Decoherence  Time  T2  Measurements . 

14.4.1  Measurement  Procedurel . 

14.4.2  Resultsl . 

4.5  Spin-Lattice  Relaxation  Time  T\  Measurements . 

14.5.1  Measurement  Procedurel . 

14.5.2  Results  and  Limitations! . 

14.6  Conclusion! . 


5  Analysis  and  Outlook 

15.1  Qverviewl . 

5.2  Sensitivity  Figure  of  Merit| . 

5.2.1  Modelling  the  Sensitivity! .  .  . 

5.2.2  Sensitivity  vs.  Irradiation  Dose 

15.3  Outlookl . 


54 

55 

55 

57 

59 

62 

64 

67 

69 

69 

70 

71 

71 

73 

74 

74 

76 

78 

78 

79 

80 

81 

81 

81 

81 

84 

86 


12 


A  Experimental  Setup  and  Calculations 


91 


A.l  Experimental  Setup . 

A.  1.1  Overview  of  the  Setup|  .  .  . 

IA.,2  Fluorescence  Cal.cul.ati.on.sl . 

A. 2.1  Power  Density! . 

IA.2.2  Saturation  Powerl . 

A. 2. 3  Imaging  Region  Calculations 


IB  NV  Fluorescence  Methods  and  Calculations! 

B.l  Charge  State  Data  Acquisition . 

IB.  1.1  Measurement,  Procedurel . 

B.l. 2  Linear  Combination  of  NV  Spectra]  . 


91 

91 

94 

95 

96 

96 

97 

97 

97 

100 


1C  Coherence  Time  Measurement  Procedurel 

1(71  Electronicsl . 

C.2  Static  Magnetic  Field  Alignment  . 

IC.2.1  Initial  Considerations! . 

C.2. 2  (100)  Diamond  Geometry  and  Initial  Field  Alignment 

1(72.3  CAW-QDMR,  Measurements! . 

C.2.4  Laser  Polarization-Assisted  Field  Alignment! . 

C.3  Resonant  Excitation  and  Field  Projection  . 

1(74  Rabi  QscillationsI . 

C.5  Polarization  Time  Optimization! . 


105 

105 

106 
106 

107 

108 
109 
112 

114 

115 


13 


14 


List  of  Figures 


1-1  The  nitrogen-vacancy  center  in  diamond 


1-2  Electronic  energy  level  structure  of  the  NY  center  in  diamond 


27 

29 


1-3  Photoluminescence  of  the  NV  center,  (left)  Diagram  depicting  the  optical  and 
vibrational  transitions  that  occur  under  optical  excitation.  Non-radiative  de- 
cay  by  vibrational  states  coupled  to  the  main  3A2  -H-3  E  NV~  transition  result 
in  a  broad  emission  pattern  from  637  -  800  nm  at  room  temperature,  (right) 
Photoluminescence  spectra  of  the  NV°  and  NV~  charge  states,  normalized  to 
I  have  the  same  areal .  30 


1-4  Currently  proposed  NV°  energy  level  structure.  The  575  nm  ZPL  transition 
is  between  a  proposed  spin  doublet  ground  state  and  an  optically  excited  spin 
doublet  excited  state.  Experimental  evidence  also  suggests  the  existence  of  a 
metastable  spin  quartet  stateT] .  32 


1-5  Two-photon  charge  state  conversion  between  the  negative  and  neutrally  charged 
NV  center.  At  wavelengths  below  477  nm  ionization  from  NV~  to  NV°  is  pre- 
dominantly  a  one-photon  process,  as  with  recombination  from  NV°  to  NV~ 
below  422  nm.l .  33 


1-6  Spin-state  dependent  dynamics  of  the  NV  center.  The  ±  1)  excited  states 
exhibit  a  higher  branching  ratio  to  decay  via  the  non-radiative  ISC  transition. 


After  decay  to  the  1E  state,  the  0)  ground  state  is  preferentially  populated, 


leading  to  spin  polarization.  The  spin-dependent  ISC  and  optical  polarization 
lead  to  spin-state  dependent  fluorescence,  where  the  emitted  light  intensity 


under  optical  excitation  maps  out  the  NV  spin  state.  .  34 


15 


1-7  Basic  detection  of  an  applied  magnetic  field  through  a  Zeeman  splitting  of 


ODMR  spectra,  (a)  The  NV  fluorescence  is  monitored  as  MW  excitation  is 


swept;  when  the  MW  is  on  resonance  with  a  transition  between  the  different 


spin  projections  of  the  3A2  state,  the  fluorescence  decreases,  (b)  An  applied 


magnetic  field  splits  the  degeneracy  between  the  ±1)  states,  giving  rise  to  two 


ODMR  dips.  The  distance  between  the  dips  is  proportional  to  the  projection 


of  the  field  along  the  NV  axis. .  35 


1-8  Modeled  fluorescence  data  that  is  typically  seen  in  magnetometry  measure- 
ments.  The  contrast  a  in  both  curves  corresponds  to  the  deviation  from  the 
maximum  fluorescence  level  that  is  caused  by  applied  MW  fields  and  the  flu¬ 


orescence  /?  corresponds  to  the  luminescence  signal  from  the  NV  ensemble. 
For  the  DC  (left)  scenario  T2*  is  a  measure  of  the  broadening  of  the  lineshape, 


while  in  the  AC  (right)  scenario  T2  characterizes  the  decay  envelope 


38 


2-1  The  effects  of  diamond  processing  on  NV  ODMR  lineshapes.  (Left)  A  diagram 
of  the  diamond  lattice  showing  a  nitrogen-vacancy  center  (circled  in  green) 
along  with  other  defects  (circled  in  red)  that  interact  with  the  NV  center  and 
broaden  its  resonance  profile,  such  as  single  vacancies,  substitutional  nitrogen 
atoms,  or  13C  atoms.  (Right)  Simulated  ODMR  curves  for  (top)  a  naturally 
abundant  98.9%  12C  diamond  and  (bottom)  a  isotopically  12C  pure  diamond, 


where  the  effects  of  broadening  from  13C  are  eliminated 


42 


2-2  Three  Element  Six  standard  grade  diamonds  used  in  this  thesis  with  the  same 
electron  dose  subjected  to  varying  annealing  temperatures.  The  change  in 
color  is  due  to  the  annealing  out  of  different  defects  at  higher  temperatures.|  50 


3-1  Example  photoluminescence  spectra  for  Diesel  (1  x  1015  e/cm2,  1250°C  an¬ 
neal)  at  10  mW  of  laser  excitation.  The  PL  signatures  from  the  NV°  and 
NV~  charge  states  are  also  shown  to  highlight  the  contribution  of  each  charge 


state,  as  well  as  the  overlap  of  the  two  spectra 


53 


16 
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3-12  Photoluminesce  spectra  of  the  diamond  Cannon  (5  x  1019  e  /cm2,  850°C 

anneal)  at  all  laser  powers.  Contrary  to  diamonds  irradiated  at  a  lower  dose, 

the  NV  contribution  to  the  PL  spectra  increases  as  a  function  of  laser  power. 

66 

4-1  (top)  Ramsey  pulse  sequence  used  to  measure  spin-dephasing  time  T2*.  (bot¬ 
tom)  Example  Ramsey  curve  from  the  diamond  Rutherford  (1  x  1016  e“/cm2, 
1250°C  anneal),  where  the  calculated  spin-dephasing  time  is  T|  =  719  ±  41.6 
ns. I .  72 


4-2  Spin-dephasing  times  T|  as  a  function  of  dose  measured  with  the  Ramsey  pulse 


sequence.  Due  to  environmental  factors  that  were  not  explicitly  controlled, 
these  values  are  to  be  seen  as  a  lower  bound  of  the  intrinsic  T2*  of  the  diamond, 
not  limited  by  external  environmental  factors,  the  dotted  line  indicate  the 


expected  13C-limited  T2  =  1  /rs  [lj 
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4-3  (top)  Pulse  sequence  used  to  measure  the  spin-decoherence  time  T2.  (bottom) 
Example  signal  from  the  diamond  Noddack  (5  x  1016  e“/cm2,  1250°C  anneal), 
where  the  calculated  spin-decoherence  time  is  T2  =  450  ±2.59  /is.  The  col- 
lapses  and  revivals  in  the  signal  are  due  to  the  Larmor  precession  of 13  C  nuclei 
in  the  surrounding  environment.! .  75 

4-4  Spin-decoherence  times  T2  for  all  diamonds  as  a  function  of  irradiation  dose.|  77 


4-5  (top)  Pulse  sequence  used  to  measure  the  spin- relaxation  time  Ti.  (bottom) 

Example  signal  from  a  spin-relaxation  pulse  sequence,  where  the  calculated 

spin-relaxation  time  is  T\  =  5.4  zb  0.0751  ms 
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4-6  Spin-lattice  relaxation  times  J\  as  a  function  of  dose.  All  diamonds  that 
had  Ti  measurements  under  high  laser  excitation  and  experienced  partial  re¬ 


polarization  during  relaxation  are  omitted 


79 


5-1  Sensitivity  figure  of  merit  FOM  as  a  function  of  long-pas  filter  wavelength  A0 


for  Rutherford  (1  x  1016  e  /cm2,  1250°C  anneal)  at  70  mW  of  excitation. 

.  .  84 

5-2  Optimum  long-pass  filter  wavelength  A0  as  a  function  of  laser  power  and  irra- 

diation  dose  for  diamonds  annealed  at  (top)  850°C  and  (bottom)  1250°C. 

.  .  85 
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5-3  Relative  sensitivity  vs.  irradiation  dose  over  four  orders  of  magnitude  of  exci- 


tation  power  for  diamonds  that  underwent  (a)  850° C  annealing  treatment  (b) 


1250°C  annealing  treatment. 
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5-4  Relative  sensitivity  vs.  irradiation  dose  at  all  anneal  temperatures  for  (left) 


100  /iW  and  (right)  150  mW  of  laser  excitation .  87 


5-5  Example  EPR  spectrum  for  Ladd-Franklin  (5  x  1018  e  /cm2,  850°C  anneal). 
While  the  signal  ideally  should  have  a  flat  baseline,  the  MW  power  causes  a 


response  from  holding  tube  that  needs  to  be  subtracted. 
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A-l  Diamond  characterization  experimental  setup  (top)  photograph  and  (bottom) 


schematic! .  92 


A-2  Diamond  mount,  (left)  Vertical  view  of  the  objective-diamond  assembly. 
Laser  light  reflects  up  to  the  objective  in  a  cage  mount  where  it  focuses  onto 


the  diamond  mount  translation  stage,  which  is  adjustable  to  ensure  that  the 
confocal  volume  is  at  the  top  of  the  diamond  near  the  MW  loop  for  coher¬ 
ence  measurements  or  near  the  bottom  of  the  diamond  for  PL  spectroscopy 
measurements,  (right)  Left:  photograph  of  the  diamond  mount.  The  loop 
antenna  is  pressed  against  the  face  of  the  diamond  and  centered  about  the 
confocal  spot  to  maximize  MW  uniformity.  Right:  diamond  mount  schematic 


with  top  and  side  views 
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A-3  Profile  of  the  laser  beam  right  before  being  reflected  off  the  dichroic.  The 


D4crX  and  D4crY  give  measures  of  the  beam  widths  in  the  X  and  Y  directions. 
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B-l  Experimental  apparatus  for  PL  spectroscopy  (a)  Spectrometer.  A  fiber  is 
mounted  to  the  entrance  slit  where  the  acquired  NV  fluorescence  is  diffracted 


and  processed  by  the  CCD  camera.  The  position  of  the  fiber  can  be  adjusted 
in  two  dimensions.  A  micrometer  allows  fine  adjustment  of  the  spectrometer 


input  slit,  (b)  Fiber  shaker  to  eliminate  modal  noise  in  the  acquired  PL  spectra.  98 
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B-2  Spectrum  modification  to  correct  for  experimental  artifacts,  (top)  Normalized 
photoluminescence  spectra  of  the  diamond  Hilbert  (5  x  1015  e/cm2,  1250°C 
anneal)  at  1  mW  optical  excitation  with  and  without  the  intensity  calibration 
turned  on.  (bottom)  Photoluminescence  spectra  of  the  diamond  Anderson 
(1.6  x  1019  e“/cm2,  1250°C  anneal)  at  300  /iW  optical  excitation  with  the 
fiber  shaker  turned  off  and  on.  Notice  how  the  NV~  zero-phonon  line  is 
I  resolvable  with  the  fiber  shaker  onl .  99 


B-3  (Left)  PL  spectra  for  Bascom  (5x  1017  e/cm2,  850°C  anneal)  at  200  mW  (top) 
and  Heaviside  (unprocessed)  at  100  /IW  (bottom).  Bascom  serves  as  the  NV° 
basis  function  due  to  its  high  fluorescence  and  negligible  NV-  content.  (Right) 
Construction  of  the  NV~  basis  function.  The  NV°  basis  function  is  scaled  to 

fit  the  NV°  ZPL  of  Heaviside’s  PL  spectrum,  which  is  then  subtracted  to  yield 
the  NV~  basis  function.  The  original  and  basis-fitted  Heaviside  spectra  are 
I  also  shown. I .  101 


B-4  PL  spectra  for  Bascom  (5  x  1017  e  /cm2,  850°C  anneal)  at  200  mW  for  405 


nm  and  532  nm  excitation^ 
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0-1  Wiring  diagram  for  coherent  control  of  the  NV  center.  PBO  controls  collection 
of  photons;  PB2  controls  the  MW  field  delivered  to  the  diamond;  and  PB3 


switches  the  laser  to  excite  and  polarize  the  NV  center 


106 


C-2  Diagram  of  the  dipole  magnet  and  (100)  diamond  geometry.  (Left)  Top  view 
of  the  magnet  and  diamond.  The  diamond  is  mounted  in  the  setup  such  that 
the  corner  is  bisected  by  the  dipole  axis.  (Right)  Side  view  of  the  diamond 
and  magnet.  The  ideal  alignment  for  this  setup  and  diamond  orientation  is 


to  have  the  field  axis  aligned  along  the  NV  class  shown 


108 
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C-3  Continuous- Wave  Optically  Detected  Magnetic  Resonance,  (top)  Pulse  se¬ 
quence  to  obtain  CAW-ODMR  data,  (bottom)  Measured  CAW-ODMR  spec¬ 
trum  of  the  diamond  Noddack  (5  x  1016  e/cm2,  1250°C  anneal).  The  signal 
mi /mo  shows  the  deviation  from  maximum  fluorescence  when  the  MW  drives 
I  an  NY  orientation  at  or  near  resonance.  This  data  was  taken  when  the  static  I 

field  was  aligned  to  one  NV  class  (outermost  resonance  pair)  and  triply  degen¬ 
erate  with  the  other  three  classes  (innermost  pair) .  The  applied  MW  power  is 
set  low  to  observe  the  hyperfine  splitting  of  each  resonance  due  to  interaction 
with  the  14N  nucleus. I .  110 


C-4  The  effect  of  laser  polarization  on  NV  class  excitation.  At  one  angle  (130° 
with  the  axes  set  by  the  third  HWP  in  the  setup),  the  green-colored  pair 
is  maximally  excited  as  the  electric  field  is  perpendicular  to  the  NV  axis. 
Rotating  the  HWP  45°  in  either  direction  maximally  excites  the  other  NV 
pair.  Selectively  exciting  one  pair  over  another  aids  in  magnetic  field  alignment. [Ill 


C-5  CAW-ODMR  spectrum  of  all  four 


0)  <r> 


1)  transitions  for  the  diamond 


Hilbert  (5  x  1015  e“/cm2,  1250°C  anneal)  for  different  half-wave  plate  angles. 
The  magnetic  field  is  mis-aligned  here  such  that  all  four  resonances  are  distin¬ 
guishable.  Inset:  maximum  deviation  from  the  reference  fluorescence  of  the 


left-most  NV  class  as  a  function  of  the  half-wave  plate  angle  .| . 112 


C-6  CAW-ODMR  spectra  of  the  (left)  0)  -B-  —  1)  and  (right)  0)  -B  1)  elec¬ 


tronic  transitions  for  the  diamond  Ramon  y  Cajal  (1.6  x  1018  e/cm2,  850°C 
anneal)  The  three  Lorentzian  peaks  arise  from  the  hyperfine  interaction  with 
the  14N  nucleus,  which  can  only  be  resolved  when  the  NV  ensemble  is  not 


power-broadened  from  the  MW  and  optical  drive.  The  slight  polarization  of 
the  peaks  closest  to  the  zero-field  splitting  of  2.87  GHz  arises  from  partial 


polarization  of  the  nuclear  spin 
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C-7  Rabi  nutation  experiment,  (top)  Pulse  sequence  used  to  measure  popula¬ 
tion  oscillations  under  a  resonantly  applied  MW  field,  (bottom)  Fluorescence 
contrast  mi/rriQ  as  a  function  of  MW  pulse  duration  for  the  diamond  Bohr 
(5  x  1017  e/cm2,  1250°C  anneal).  The  coherent  population  oscillations  yield 


a  sinusoidal  huorescence  signal  with  a  stretched  exponential  decay  envelope.|  116 


C-8  Optical  pumping  pulse  length  determination  (top)  Pulse  sequence  used  to 


measure  the  re-pump  time  necessary  to  optically  polarize  the  NY  into  the  0) 


state,  (bottom)  Example  signal  from  a  pump  measurement  for  the  diamond 


Lorentz  (5  x  1018  e"/cm2,  1250°C  anneal).  After  optical  pumping,  a  n  pulse 
hips  the  spin  population  and  a  repump  laser  pulse  is  applied.  After  a  time  r, 
the  fixed-length  gate  pulse  is  turned  on  to  measure  the  huorescence  m\.  The 


sequence  is  repeated  without  the  7r  pulse  to  give  the  signal  m0.  After  sweeping 
the  start  time  r,  the  time  after  which  both  signals  achieve  their  steady  state 
values  is  chosen  as  the  minimum  duration  of  the  repump  pulse| . 117 
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Chapter  1 


Magnetometry  and  the 
Nitrogen- Vacancy  Center 

1.1  Introduction 

The  measurement  of  increasingly  precise  magnetic  fields  has  served  to  give  civilizations 
the  strategic  edge  with  regards  to  military  warfare,  navigation/exploration,  and  scientific 
investigation.  The  earliest  recorded  use  of  a  magnetometer  was  by  the  Chinese  emperor 
Hoang-Ti  (~2700  B.C.E.),  when  he  navigated  his  army  through  fog-shrouded  lowlands  using 
a  magnetic  rock,  called  loadstone,  to  capture  enemy  forces  [2].  This  same  loadstone  served 
as  the  magnetization  source  to  “feed”  iron  needles  in  compasses  as  16th  century  explorers 
mapped  the  new  world  in  the  Age  of  Exploration.  In  the  19th  and  20th  centuries,  advances  in 
classical  and  quantum  physics  ushered  in  new  generations  of  highly-sensitive  magnetometers 
that  continue  to  push  the  limits  of  human  exploration.  Modern  applications  today  include 
mapping  the  magnetic  field  of  the  brain  [3J  and  assisting  archaeologists  in  magnetically 
mapping  temples  of  ancient  civilizations  [4J.  The  push  for  higher  sensitivity  in  magnetic 
field  detection  and  imaging  promises  to  unearth  new  applications  in  the  life  and  material 
sciences  in  addition  to  increasing  our  understanding  of  mesoscopic  physics. 

Certain  limits  exists  for  today’s  magnetometers,  however.  While  atomic  magnetometers 
currently  possess  the  highest  sensitivity  in  magnetometry  |5J,  these  systems  typically  measure 
only  the  magnitude  of  the  magnetic  field;  i.e.  they  are  scalar  magnetometers.  To  measure 
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both  field  magnitude  and  direction,  vector  magnetometers  (such  as  SQUIDS  and  fluxgates) 
rely  on  mechanical  axes  of  the  system,  which  can  be  limited  by  mechanical  and  electrical 
drifts.  Recently,  the  nitrogen-vacancy  (NV)  center  in  diamond  has  emerged  as  a  favorable 
candidate  for  high-sensitivity  scalar  and  vector  quantum  metrology  mm-  In  the  last  five 
years,  ensembles  of  NV  centers  have  mapped  the  //T-fields  in  /rm-size  meteorite  grains  |8J  and 
noninvasive  mapped  the  action  potentials  of  single  neurons  [9j.  In  addition,  NV  ensembles 
in  diamond  are  inherent  vector  magnetometers  with  the  advantage  of  measurement  axes  tied 
to  the  unchanging  crystal  axes  of  diamond. 

An  active  area  of  investigation  is  the  improvement  of  the  NV  spin-state  dependent 
fluorescence-based  measurement  of  magnetic  fields.  This  thesis  presents  a  systematic  study 
of  the  effects  of  diamond  processing  to  understand  the  effects  of  irradiation  and  annealing  on 
NV  formation  and  the  spin  properties  in  commercially  available  bulk  chemical  vapor  deposi¬ 
tion  (CVD)  diamond  for  ensemble  magnetometry.  Nitrogen- vacancy-center  charge  state  and 
spin  decoherence  times  are  measured  for  diamonds  irradiated  with  1  ME  electrons  at  doses 
of  1  x  1015  -  5  x  1019  e  /cm 2  and  subsequently  thermally  annealed  at  temperatures  of  850°  C 
and  1250°C.  The  study  provides  an  optimal  range  for  diamond  processing  and  shows  the 
quenching  of  the  NV  center  at  high  dosage  from  the  creation  of  additional  vacancy-related 
defects. 


1.2  Nitrogen- Vacancy  Center  Physics 

Among  the  hundreds  of  color  centers  in  diamond  nm  is  the  nitrogen-vacancy  center,  which 
is  composed  of  a  substitutional  nitrogen  atom  adjacent  to  a  vacancy  in  the  diamond  lattice. 

shows  this  defect  along  the  [111]  crystallographic  axis  in  diamond,  where  it  exhibits 
trigonal  symmetry  about  this  axis.  The  NV  center  is  primarily  found  in  one  of  two  charge 
states:  the  negatively  charged  NV-  state  and  the  neutrally  charged  NV°  state.  The  NV- 
charge  state  is  the  focus  of  this  thesis  due  to  its  favorable  spin  and  optical  properties  for 
high-sensitivity  magnetometry.  The  NV°  charge  state  exhibits  markedly  different  properties 
than  the  NV-  center  and  is  not  suited  for  magnetometry;  as  will  be  explained  below,  it 
can  actually  hinder  the  operation  of  an  NV  magnetometer.  In  this  thesis,  unless  explicitly 


Figure  1-1 
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Figure  1-1:  The  nitrogen- vacancy  center  in  diamond. 

stated,  the  term  NV  center  will  refer  to  the  negative  charge  state  NV-. 

The  first  observation  of  the  NV  color  center  came  in  1965  when  du  Preez  m  observed 
the  absorption  spectrum  of  a  high-nitrogen  content  diamond  subject  to  0.78  ME  electron 
irradiation  and  thermal  annealing  at  temperatures  exceeding  600°C.  As  will  be  discussed  later 
in  this  thesis,  irradiation  from  highly  energetic  particles  creates  vacancies  in  the  diamond 
lattice,  while  thermal  annealing  causes  the  vacancies  to  diffuse  throughout  the  diamond.  The 
vacancies  are  trapped  by  the  nitrogen  to  form  the  NV  center.  Systematically  increasing  the 
annealing  temperature  and  measuring  the  absorption  spectra  after  each  anneal,  after  600°C 
du  Preez  noticed  the  optical  bands  in  the  absorption  spectra  corresponding  to  the  negative 
vacancy  (ND1)  and  neutral  vacancy  (GR1)  decreased  as  a  new  band,  the  “6400  A”  (640  nm) 
optical  line  increased.  This  of  course  is  the  well-known  637  nm  optical  transition  of  the  NV- 
center.  These  considerations  led  du  Preez  to  speculate  the  band  belonged  to  a  substitutional 
nitrogen-vacancy  pair  in  the  diamond  lattice. 

In  the  1970s  the  physical  and  electronic  structure  of  the  “640  nm”  band  was  unraveled. 
Photoluminescence  measurements  taken  by  Clark  and  Norris  [12J,  in  addition  to  the  uniax¬ 
ial  stress  measurements  of  Davies  and  Hamer  [13],  established  the  trigonal  symmetry  of  the 
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defect  and  elucidated  its  electronic  structure.  Electron  paramagnetic  resonance  (EPR)  mea¬ 
surements  by  Loubser  and  van  Wyk  led  to  the  assignment  of  the  negative  charge  state  and  a 
molecular  orbital  model  to  model  the  electronic  structure.  In  this  model,  the  dangling  bonds 
of  the  carbon  and  nitrogen  atoms  surrounding  the  vacancy  provide  five  of  the  six  electrons 
attributed  to  NV-,  while  the  sixth  is  assumed  to  come  from  a  nearby  charge  donor.  The 
current  explosion  in  the  interest  of  the  NV  center,  however,  can  be  traced  back  to  1997,  when 
Gruber  et.  al.  used  the  evolving  techniques  of  confocal  microscopy  to  optically  address  and 
observe  single  NV  centers  in  diamond  M-  It  was  quickly  discovered  afterward  that  the  NV 
center  possessed  long  coherence  times  and  that  its  spin  state  could  be  readout  with  conven¬ 
tional  optically-detected  magnetic  resonance  techniques.  Among  the  many  applications  the 
NV  is  suited  for,  these  properties  endow  the  NV  center  as  a  highly-sensitive  magnetometer. 

This  section  will  briefly  overview  the  energy  level  structure  of  the  NV  center,  to  include 
the  neutrally-charged  NV°.  The  spin  polarization  and  spin-state  dependent  fluorescence 
emission,  two  crucial  ingredients  for  practical  use  of  the  NV  center,  will  then  be  explained. 
A  more  in-depth  description  of  the  NV  centers  can  be  found  in  [T5J. 


1.2.1  Electronic  Structure 


Figure  1-2  shows  a  simplified  diagram  of  the  NV  center’s  electronic  structure.  We  will  only 


examine  the  room  temperature  energy  structure  in  this  thesis.  A  full  group  theoretical 
treatment  of  the  NV  center  [16J  reveals  that  a  set  of  molecular  orbitals  form  the  basis  of 
the  NV  energy  structure.  These  are  deep  on  the  diamond  band-gap  and  are  filled  with  the 
six  electrons  of  the  NV  center.  Instead  of  thinking  in  terms  of  six  electrons,  however,  it  is 
simpler  to  think  of  the  energy  structure  in  terms  of  two  missing  electrons  (“holes”)  that  are 
needed  to  fill  the  molecular  orbitals.  These  two  holes  give  rise  to  a  spin-triplet  ground  state 
and  excited  state.  Here  we  take  the  spin  projection  axis  to  be  the  NV  axis.  The  electronic 
ground  state  is  a 1 *  3A2  spin  triplet^  [10J  with  a  zero  field  splitting  of  Dgs  =  2.87  GHz  between 
the  |0)  ( ms  =  0)  and  |  ±  1)  ( ms  =  ±1)  spin  sublevels  fTTTj .  This  splitting  is  caused  by  the 


1The  labels  A1,  A2,  and  E  describe  the  three  irreducible  representations  of  the  C:>v  symmetry  point  group 

of  the  NV  center.  Each  of  these  has  a  certain  symmetry  and  constrain  the  types  of  dip  ole- allowed  transitions 

that  can  take  place. 
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Figure  1-2:  Electronic  energy  level  structure  of  the  NV  center  in  diamond. 

spin-spin  interaction  of  the  electrons.  The  excited  state  3  E  is  also  a  spin  triplet  state  with 
a  zero  held  splitting  of  Des  =  1.42  GHz  D3-  The  |0)  states  of  the  ground  and  excited  spin 
triplet  states  are  separated  by  a  1.945  eV  (637  nm)  zero-phonon  line  (ZPL).  Application  an 
external  magnetic  held  B  causes  a  Zeeman  shift  between  the  |±1)  levels,  leading  to  an  energy 
splitting  of  2'yiwB\\.  Here  yyvv  is  the  NV  gyromagnetic  ratio  (2.8  MHz/G  [15j)  and  B\\  is 
the  component  of  the  magnetic  held  along  the  NV  axis.  The  nitrogen  nuclear  spin  causes 
an  additional  hyperhne  splitting  of  each  electronic  spin  level  into  two  (for  15N,  I  =  1/2)  or 
three  (14N,  1  =  1)  substates.  Finally,  in  addition  to  the  spin  triplets  are  two  intermediate 
singlet  states  3A\  and  1E  separated  by  a  1.190  eV  (1042  nm)  ZPL  [18j. 

1.2.2  Photoluminescence  Spectra 

Both  the  NV°  and  NV-  have  resonant  transitions  at  their  respective  zero-phonon  lines. 
However,  the  optimal  excitation  wavelength  of  the  NV-  center  (to  minimize  population  in 
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Figure  1-3:  Photoluminescence  of  the  NV  center,  (left)  Diagram  depicting  the  optical  and  vi¬ 
brational  transitions  that  occur  under  optical  excitation.  Non-radiative  decay  by  vibrational 
states  coupled  to  the  main  3A2  -f-)-3  E  NV-  transition  result  in  a  broad  emission  pattern 
from  637  -  800  nm  at  room  temperature,  (right)  Photoluminescence  spectra  of  the  NV°  and 
NV-  charge  states,  normalized  to  have  the  same  area. 


NV°  under  optical  illumination)  is  510  -  540  nm  [19j.  The  use  of  off-resonant  higher  energy 
light  excites  the  NV  into  the  continuum  of  vibrational  states  induced  by  the  vibrational 
degrees  of  freedom  of  the  electrons  and  nuclei.  Figure  1-3  (left)|  shows  the  effect  of  this 
vibronic  coupling.  The  vibrational  states  can  approximately  be  regarded  as  states  in  a 
harmonic  well,  where  the  lowest  state  corresponds  to  zero  phonons,  the  first  state  having 
coupling  to  one  phonon,  etc.  Off-resonant  light  excites  the  NV  into  a  higher-energy  phonon 
state.  Phonon  relaxation  quickly  brings  the  NV  spin  back  to  one  of  the  electronic  states, 
where  it  can  then  spontaneously  decay  to  the  electronic  ground  state  through  a  zero-phonon 
transition  (575  nm  for  NV°,  637  nm  for  NV-)  or  to  a  higher-energy  phonon  state,  resulting 
in  a  range  of  longer  wavelengths  called  the  phonon-sideband  (PSB).  The  PSB  is  the  primary 
decay  channel  for  the  NV  center;  for  the  NV-  center  the  emission  into  the  ZPL  is  only  a  few 
percent  m- 


depicts  the  photoluminescence  (PL)  spectrum  of  NV°  and  NV-,  which 
is  the  intensity  of  light  emitted  as  a  function  of  wavelength.  The  narrow  peaks  at  575  nm 
and  637  nm  correspond  to  the  zero-phonon  lines,  whereas  the  broad  emission  band  at  higher 


Figure  1-3  (right) 
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wavelengths  corresponds  to  the  PSB.  Note  the  overlap  of  the  NV°  PL  spectrum  with  that 
of  NV-.  We  will  see  that  this  adds  background  fluorescence  to  NV-  state  readout,  which  is 
a  problem  in  magnetometry  measurements. 


1.2.3  The  Neutrally-Charged  NV  Center 


As  mentioned  above,  the  other  common  charge  state  of  the  NV  center  is  the  neutrally- 
charged  NV°  center.  Initially  observed  in  cathodoluminescence  spectra  in  1960  m,  this 
center  was  detected  through  optical  absorption  and  emission  measurements  in  1979  by  Davies 
m ■  Conclusive  assignment  of  the  observed  2.156  eV  optical  band  to  NV°  occurred  in  1996 
after  Mita  [23]  noticed  the  correlated  changes  in  NV-  and  NV°  fluorescence  under  increased 
neutron  irradiation.  At  first  it  was  proposed  that  changes  in  the  Fermi  level  from  diamond 
impurities  were  responsible  for  the  charge  state  of  the  NV  center,  but  Collins  [24J  argued 
that  this  would  not  permit  the  co-existence  of  both  charge  states.  Instead,  the  microscopic 
distribution  of  proximal  charge  donors  plays  the  dominant  role  in  influencing  the  NV’s  charge 
state. 

Since  it  has  one  less  electron  than  NV-,  the  electronic  energy-level  structure  of  NV° 
differs  from  that  of  NV-.  |Figure  P4  shows  the  currently  accepted  energy  level  structure. 
Note  that  the  NV°  energy-level  structure  of  NV°  is  less  studied  than  NV-.  The  proposed 
NV°  states  are  a  spin  doublet  ground  and  optically  excited  state,  in  addition  to  a  spin-quartet 
metastable  state  [25j.  Since  the  NV°  does  not  exhibit  the  same  spin-dependent  fluorescence 
as  NV-  under  MW  excitation,  the  overlapping  PL  of  NV°  adds  background  fluorescence  that 
degrades  the  signal-to-noise  of  NV-  fluorescence  measurements. 

It  is  thus  desirable  for  magnetometry  purposes  to  eliminate  the  influence  of  NV°  or  to  use 
it  for  readout.  In  general,  elimination  of  NV°  fluorescence  is  difficult  as  several  phenomena 
contribute  to  charge  conversion  processes  in  addition  to  the  native  NV°  population  in  the 
lattice.  With  regards  to  the  diamond  composition,  if  the  vacancy  concentration  is  comparable 
to  the  nitrogen  concentration  in  the  diamond  lattice,  NV°  will  preferentially  form  over  NV-. 
As  explained  above,  this  is  due  to  influence  the  local  environment  of  a  nitrogen- vacancy  has 
in  determining  the  charge  state.  In  particular,  if  the  single  substitutional  nitrogen  atoms  are 
increasingly  paired  to  vacancies,  the  amount  available  to  donate  an  electron  to  an  existing 
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Figure  1-4:  Currently  proposed  NV°  energy  level  structure.  The  575  nm  ZPL  transition  is 
between  a  proposed  spin  doublet  ground  state  and  an  optically  excited  spin  doublet  excited 
state.  Experimental  evidence  also  suggests  the  existence  of  a  metastable  spin  quartet  state. 


NV°  center  to  create  NV  decreases. 

Even  if  an  NY-  center  exists  in  diamond  after  growth  and  processing,  under  optical 


illumination  the  NV-  center  can  be  ionized  to  form  the  NV°  center.  Figure  1-5  shows  the 
two-photon  ionization  and  recombination  charge  state  dynamics  [19j  between  the  NV-  and 
NV°.  Starting  in  the  NV-  center,  under  optical  excitation  (1)  an  electron  in  the  NV-  is 
promoted  to  its  excited  state  by  one  photon.  (2)  From  there,  electron  is  then  excited  to  the 
conduction  band  of  the  host  diamond  by  a  second  photon,  and  the  resulting  NV°  center  is 
left  in  either  its  excited  state  (which  rapidly  decays)  or  ground  state.  In  the  recombination 
process,  (3)  a  single  photon  first  excites  an  electron  from  the  NV°  ground  state,  and  (4) 
a  second  photon  captures  an  electron  from  the  valence  band,  leaving  the  NV  center  in  the 
ground  state  of  NV-.  Under  common  operating  conditions  with  a  low  intensity  green  laser, 
the  steady-state  occupation  probability  for  the  NV  to  be  in  the  negative  charge  state  is 
~  70%  [26J.  Note  that  at  high  enough  photon  energies  NV-  can  be  directly  ionized  to  NV° 
with  a  single  photon,  which  increases  the  probability  to  be  in  NV°.  At  still  higher  energies 
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Figure  1-5:  Two-photon  charge  state  conversion  between  the  negative  and  neutrally  charged 
NV  center.  At  wavelengths  below  477  nm  ionization  from  NV-  to  NV°  is  predominantly  a 
one-photon  process,  as  with  recombination  from  NV°  to  NY-  below  422  nm. 


the  recombination  process  can  be  accomplished  with  a  single  photon,  thus  increasing  the 
NV-  occupation  probability  127J. 


1.2.4  Spin-State  Dependent  Fluorescence  and  ODMR 

We  conclude  this  section  by  discussing  two  physical  mechanisms  that  most  applications  of 
the  NV-  center  are  critically  dependent  on.  When  the  NV-  center  is  optically  excited, 
there  are  two  primary  ways  for  the  excited  state  to  decay,  depending  on  the  initial  electronic 
ground  state.  While  the  main  decay  path  is  to  spontaneously  decay  through  the  zero-phonon 
transition  at  637  nm  or  into  the  phonon-sideband  between  640  -  800  nm,  there  exists  a  non- 
radiative  decay  pathway  via  the  intersystem  crossing  (ISC)  to  the  XA\  singlet  state.  The 
weak  coupling  between  the  3E  and  lA\  states  arises  from  spin-orbit  coupling;  this  also  weakly 
couples  the  1E  and  3A2  manifolds.  If  the  NV-  spin  decays  through  the  ISC,  it  quickly  decays 
to  1E  by  emission  of  an  infrared  photon,  which  then  decays  to  the  ground  state  3A2. 

Two  properties  are  of  interest  here:  first,  the  optical  cycling  between  the  ground  and 
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Figure  1-6:  Spin-state  dependent  dynamics  of  the  NV-  center.  The  |  ±  1)  excited  states 
exhibit  a  higher  branching  ratio  to  decay  via  the  non-radiative  ISC  transition.  After  decay 
to  the  1E  state,  the  |0)  ground  state  is  preferentially  populated,  leading  to  spin  polarization. 
The  spin-dependent  ISC  and  optical  polarization  lead  to  spin-state  dependent  fluorescence, 
where  the  emitted  light  intensity  under  optical  excitation  maps  out  the  NV-  spin  state. 


excited  states  is  largely  spin  conserving,  while  the  ISC  is  mostly  spin-dependent.  The  ISC 
decay  occurs  more  frequently  from  the  |  ±  1)  excited  state  sublevels.  The  NV  system  cannot 
optically  cycle  while  in  the  lE  singlet  state  (with  an  average  lifetime  of  about  250  ns  |EB]), 
hence  it  is  dark  during  this  time  period.  Since  the  |0)  state  decays  mainly  through  optical 
cycling  (average  excited  state  lifetime  of  12  ns  |29|),  it  is  deemed  a  “bright  state”  and  has 
a  higher  average  fluorescence.  The  second  property  of  interest  stems  from  the  preferentially 
decay  from  lE  to  the  3A2  |0)  sublevel.  Upon  continuous  optical  excitation,  the  state  is 
polarized  into  the  non-Boltzmann  steady  state  spin  alignment  along  the  |0)  axis  [15j.  Both 
of  these  mechanisms  (ground  state  spin  polarization  and  spin-state-dependent  fluorescence) 
are  used  to  optically  detect  the  NV  ground  spin  state:  if  the  NV  center  is  optically  excited, 
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Figure  1-7:  Basic  detection  of  an  applied  magnetic  field  through  a  Zeeman  splitting  of  ODMR 
spectra,  (a)  The  NV  fluorescence  is  monitored  as  MW  excitation  is  swept;  when  the  MW  is 
on  resonance  with  a  transition  between  the  different  spin  projections  of  the  3A2  state,  the 
fluorescence  decreases,  (b)  An  applied  magnetic  field  splits  the  degeneracy  between  the  |  ±1) 
states,  giving  rise  to  two  ODMR  dips.  The  distance  between  the  dips  is  proportional  to  the 
projection  of  the  field  along  the  NV  axis. 


higher  fluorescence  recorded  indicates  the  spin  state  was  initially  zero.  Under  the  same 


measurement  time,  lower  fluorescence  means  the  state  was  in  |  ±  1).  Figure  1-6  shows  the 
different  transitions  involved  in  this  spin-state  dependent  fluorescence. 

The  experimental  procedure  of  optically  detected  magnetic  resonance  (ODMR)  takes 
advantage  of  the  spin-state-dependent  fluorescence  dependence  to  read  out  the  transition 


frequencies  between  the  different  spin-sublevels  of  the  NV  ground  state  (see  Appendix  C  for 
a  more  in-depth  treatment).  This  can  be  achieved  under  continuous  (CAW)  laser  excitation 
or  by  a  pulsed  technique  [30J;  here  we  will  only  discuss  the  CAW  technique.  The  NV  center 
is  continuously  excited  (usually  by  green  light),  which  causes  the  ground  state-spin  to  be 
polarized  to  |0).  This  results  in  maximum  fluorescence  due  to  the  lower  branching  ratio  to 
the  non-radiative  ISC  decay  path,  which  is  recorded.  Simultaneously,  the  frequency  of  an 
applied  MW  field  is  swept  through  a  range  of  frequencies.  When  the  MW  field  is  not  on 
resonance  between  the  |0)  <->■  |  ±  1)  ground  state  transition,  the  fluorescence  intensity  is  a 
maximum  as  the  NV  is  polarized  in  the  |0)  state.  If  the  MW  is  on  an  NV  resonance,  it  drives 
the  NV  population  which  cycles  between  the  |0)  and  |  ±  1)  states.  This  causes  a  reduction  in 
the  fluorescence  intensity  due  to  the  increased  probability  of  these  states  to  decay  using  the 
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non-radiative  ISC.  Thus,  the  recorded  fluorescence  signal  exhibits  dips  at  the  NV  resonance 
frequencies. 


In  the  absence  of  a  static  magnetic  field,  the  only  transition  detected  in  ODMR  is  the 
degenerate  |0)  -f-»  |  ±  1)  electronic  resonance  at  the  D  =  2.87  GHz  zero-field  splitting.  If  a 
static  field  is  applied  with  a  projection  B\\  along  the  NV  axis,  then  to  zeroth  order  the  |  ±  1} 
states  are  split  by  2'-/nvB\\,  where  qjw  is  the  NV  gyromagnetic  ratio  (2.8  MHz/G).  ODMR 


then  records  resonances  at  E  =  D±^nvB\\.  Figure  1-7  shows  a  simulated  fluorescence  signal 
from  an  ODMR  measurement. 


1.3  NV  Magnet ometry 

The  ODMR  technique  described  above  is  the  essence  of  the  NV  system  as  a  magnetometer, 
namely  the  detection  of  Zeeman  shifts  through  spin-state  dependent  fluorescence.  A  change 
in  the  magnetic  field  leads  to  a  change  in  the  ODMR  spectrum.  For  a  single  NV  center 
the  field  projection  along  the  NV  axis  creates  a  pair  of  Zeeman-split  energy  levels;  with 
an  NV  ensemble  as  many  as  eight  resonances  can  be  observed  due  to  the  four  different 
crystallographic  axes  in  diamond.  This  is  the  basis  of  vector  magnetometery.  The  four 
different  projections  can  be  used  to  over-determine  the  three-component  vector  field,  where 
the  axes  are  stable  to  those  of  diamond. 

In  the  magnetometry  measurement  schemes  the  central  feature  is  the  measurement  of  a 
population  difference  of  the  NV-  electronic  spin  between  |0)  and  |  ±  1).  Two  signals  are  read 
out  sequentially  in  order  to  calculate  this  population  difference.  After  optical  pumping  to 
the  |0)  state,  a  NV  fluorescence  signal  mi  is  acquired  during/after  a  prescribed  laser/MW 
excitation  protocol.  The  NV  center  is  then  re-polarized,  or  repumped,  by  the  laser  and  a 
different  MW  excitation  is  executed  to  record  a  fluorescence  signal  m0.  The  contrast  a  of 
the  measurement  is  defined  as  (m0  —  mi)/(m0  +  mi)  and  is  a  measure  of  the  population 
difference  between  the  signals  mo  and  mi.  The  number  of  collected  photons  /3  in  a  single 
measurement  is  defined  as  /3  =  (mo  +  mi)tgate/2 ,  where  tgate  is  the  fluorescence  collection 
time.  Note  that  the  measured  fluorescence  signals  m0  and  m-i  contain  fluorescence  from  NV° 
due  to  the  overlapping  PL  spectra  of  NV°  and  NV-.  However,  since  the  NV°  fluorescence 
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does  not  change  under  the  same  MW  excitation  that  changes  the  NV-  fluorescence,  the 
subtraction  of  the  two  signals  in  a  still  measures  the  NV-  population  difference,  albeit  with 
a  lower  value  due  to  the  contribution  of  the  NV°  fluorescence  in  the  denominator  of  a.  This 
compensates  for  any  improvement  in  (3  NV°  may  give. 


1.3.1  Magnetic  Sensitivity 

The  magnetic  sensitivity  in  units  of  T/  v/Hz  is  77  =  sbVt  p],  where  SB  is  the  minimum 
detectable  magnetic  field  above  the  noise  floor  (SNR  =  1)  and  the  time  T  is  the  measurement 
time.  For  NV  centers,  the  fundamental  quantum  limit  to  sensitivity  is  associated  with  spin 
projection  noise  [31j.  However,  as  most  NV  magnetometry  schemes  read  out  magnetic  fields 
through  population  differences  by  fluorescence  detection,  the  photon-shot  noise  becomes  the 
dominating  limitation  to  the  sensitivity  of  an  NV  magnetometer.  In  the  photon-shot  noise 
limited  regime,  the  magnetic  field  sensitivity  rj  takes  the  form  [32J 


7]  oc 
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VFT1 
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where  T  is  the  measurement  time,  a  is  the  fluorescence  contrast,  the  f3  is  the  average  number 
of  collected  photons.  Fundamental  constants  and  other  pre-factors  that  depend  on  the 
specific  magnetometry  protocol  are  neglected  here  for  simplicity. 

For  most  DC  magnetometry  schemes  the  optimal  sensitivity  is  obtained  when  T  is  chosen 
to  be  on  the  order  of  the  inhomogeneous  spin  dephasing  time  T2*  [30J,  while  for  most  AC 
schemes  it  is  optimized  when  chosen  to  be  on  the  order  of  the  spin  decoherence  time  T2  [7]. 
While  longer  times  in  principle  would  give  a  better  sensitivity,  when  T  is  longer  than  the 
dephasing/decohrence  time  of  the  NV  ensemble  the  contrast  of  the  sensitivity  is  degraded. 


Figure  D8|  shows  simulated  fluorescence  data  that  is  typically  measured  in  DC  and  AC 
magnetometry  experiments  and  what  each  of  the  variables  corresponds  to. 


1.3.2  Optimizing  Sensitivity 


Equation  1.1  highlights  the  three  most  important  variables  to  optimize  for  a  given  NV 


ensemble.  Physically,  a  larger  contrast  a  and  fluorescence  (3  give  a  higher  signal-to-noise  in  a 
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Figure  1-8:  Modeled  fluorescence  data  that  is  typically  seen  in  magnetometry  measurements. 
The  contrast  a  in  both  curves  corresponds  to  the  deviation  from  the  maximum  fluorescence 
level  that  is  caused  by  applied  MW  fields  and  the  fluorescence  (3  corresponds  to  the  lumi¬ 
nescence  signal  from  the  NV  ensemble.  For  the  DC  (left)  scenario  T2*  is  a  measure  of  the 
broadening  of  the  lineshape,  while  in  the  AC  (right)  scenario  T2  characterizes  the  decay 
envelope. 


magnetometry  measurement  and  can  be  improved  by  increasing  the  NV  density  in  a  diamond. 
Similarly,  a  longer  coherence  time  T2  or  T|  enables  a  longer  signal  acquisition  time  before  the 
contrast  is  degraded  by  the  surrounding  environment,  and  can  be  improved  by  having  a  high- 
purity  diamond  free  of  paramagnetic  defects  besides  the  NV  center.  Optimizing  the  three 
parameters  is  not  a  straightforward  task  as  they  are  not  independent.  For  example,  in  CAW- 
ODMR,  while  one  can  increase  the  microwave  drive  power  to  the  NV  to  increase  the  contrast 
a,  the  power  broadening  increases  the  linewidth  of  the  resonance.  Since  the  linewidth  is 
inversely  proportional  to  T2  130J,  this  compromises  sensitivity.  Similarly,  in  the  continuous 
excitation  scheme  one  can  decrease  the  laser  power  broadening  to  narrow  the  linewidth  at 
the  cost  of  decreasing  the  number  of  photons  collected.  The  same  interdependence  arises 
in  diamond  processing:  while  one  can  improve  nitrogen  to  NV  conversion  by  irradiation  of 
energetic  particles  to  increase  /3,  creation  of  NV°  and  other  defects  from  crystal  damage  can 
accompany  the  irradiation,  thus  decreasing  the  coherence  time  T2  /T2  and  contrast  a. 

We  can  now  state  the  goal  of  this  thesis:  given  the  NV  center’s  potential  for  high- 
sensitivity  magnetometry,  we  wish  to  optimize  the  sensitivity  in  order  to  detect  smaller 
fields  with  high  precision.  While  different  experimental  protocols  can  be  applied  to  achieve 
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this,  in  this  thesis  we  focus  on  the  host  diamond  itself.  That  is,  we  look  to  the  influence  of  the 
growth  and  processing  of  the  diamond  lattice  and  the  effects  it  has  on  the  sensitivity  of  an 
NV  ensemble.  We  will  characterize  a  set  of  diamonds  with  varied  processing  conditions  with 
a  set  of  measurements  to  measure  the  contrast  a,  collected  photons  fi.  and  the  coherence 
times  T2.  We  then  look  at  the  combined  effect  of  these  parameters  in  the  sensitivity  formula 
to  determine  the  processing  conditions  yielding  the  optimal  sensitivity. 

This  thesis  is  organized  as  follows:  in  |Chapter~2[  we  will  review  the  different  diamond 
processing  techniques  and  introduce  the  diamonds  used  for  this  study.  |Chapter~3l  and  [Chap¬ 


ter  4|  will  describe  the  set  of  experiments  used  to  characterize  a,  (5,  and  T  for  each  diamond. 


Finally,  all  the  experimental  results  will  be  collectively  analyzed  in  Chapter~5|  to  determine 
the  optimal  processing  parameters  for  the  specific  diamonds  used. 
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Chapter  2 


Diamond  Processing 

2.1  Overview 

The  previous  chapter  introduced  the  NV  center  as  a  unique  and  powerful  candidate  for 
scalar/vector  magnetometry  and  magnetic  field  imaging.  To  detect  sub-picoTesla  magnetic 
fields  and  to  image  sub-millimeter  magnetic  regions,  the  NV  magnetometer  must  be  tailored 
to  optimize  the  magnetic  field  sensitivity  r /.  This  sensitivity  optimization  can  be  realized, 
for  example,  by  employing  dynamical  decoupling  pulse  sequences  to  isolate  the  NV  center 
from  the  environment  [33] ,  optimizing  the  geometry  of  the  diamond  to  enhance  fluorescence 
collection  I34j,  or  by  post-diamond  growth  processing  to  enhance  the  NV  concentration  and 
purity  of  the  surrounding  environment  [35].  In  this  thesis,  we  focus  on  the  latter  approach 
and  examine  the  effects  of  diamond  processing  on  NV  magnetometry  properties.  Specifically, 
we  examine  the  effects  of  varying  electron  irradiation  doses  and  different  anneal  tempera¬ 
tures  on  a  set  of  initially  identical  diamonds  to  find  a  processing  regime  suited  to  ensemble 
magnetometry. 

Diamond  processing  in  this  context  refers  to  modifying  the  diamond  properties  during 
and  after  the  growth  process  to  produce  desired  effects.  Commercially,  these  techniques 
have  been  historically  developed  in  the  gemology  industry  where  diamonds  (among  other 
stones)  are  irradiated  with  high  energy  particles  to  change  their  color  [36] .  For  NV  diamond 
magnetometry,  the  desired  effect  is  to  increase  the  nitrogen  to  NV  conversion  efficiency 
while  also  keeping  the  environment  free  of  other  defects  that  can  decohere  the  NV  ensemble. 
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Figure  2-1:  The  effects  of  diamond  processing  on  NV  ODMR  lineshapes.  (Left)  A  diagram 
of  the  diamond  lattice  showing  a  nitrogen- vacancy  center  (circled  in  green)  along  with  other 
defects  (circled  in  red)  that  interact  with  the  NV  center  and  broaden  its  resonance  profile, 
such  as  single  vacancies,  substitutional  nitrogen  atoms,  or  13C  atoms.  (Right)  Simulated 
ODMR  curves  for  (top)  a  naturally  abundant  98.9%  12C  diamond  and  (bottom)  a  isotopically 
12C  pure  diamond,  where  the  effects  of  broadening  from  13C  are  eliminated. 


Figure  2-l]shows  the  simulated  effects  of  one  example  of  diamond  processing.  A  diamond  with 
defects  such  as  isolated  vacancies,  isolated  single  substitutional  nitrogen,  and  13C  interact 
with  the  NV  center  to  broaden  the  measured  ODMR  lineshape  |37j.  Growing  isotopically 
pure  12C  diamond  [38J  narrows  the  lineshape  by  eliminating  dipolar  coupling  with  13C,  and 
thus  increases  the  coherence  time. 


In  this  chapter  we  discuss  the  role  of  diamond  processing  in  creating  an  NV  ensemble  in 
the  diamond  lattice.  First,  two  main  methods  of  synthetic  diamond  growth  are  discussed. 
Then,  we  focus  on  the  incorporation  of  individual  nitrogen  and  vacancies  in  the  lattice,  along 
with  how  to  couple  them  to  form  an  NV  center.  Finally,  we  introduce  the  diamonds  used  in 
this  thesis,  along  with  the  processing  methods  used. 
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2.2  Diamond  Processing 


There  are  several  factors  to  consider  in  processing  a  diamond  for  NV  magnetometry:  how  to 
grow  the  diamond,  controlling  the  amount  of  nitrogen  and  vacancies  incorporated  into  the 
diamond  lattice,  and  how  to  form  the  NV  defect  centers  while  minimizing  the  presence  of 
other  unwanted  defects.  Each  has  a  significant  effect  on  the  resulting  properties  and  concen¬ 
tration  of  NV  centers  in  the  diamond  (both  charge  states)  and  the  surrounding  environment. 
This  section  will  go  over  each  of  these  factors  and  discuss  why  we  chose  specific  diamond 
processing  recipes  for  our  study. 

2.2.1  Synthetic  Diamond  Generation 

While  it  is  possible  to  perform  NV  magnetometry  with  natural  diamonds,  the  reproducibility 
and  control  of  defects  that  synthetic  diamond  affords  makes  it  the  preferred  diamond  type. 
Since  graphite  is  the  thermodynamically  stable  form  of  solid  carbon  at  room  temperature 
and  pressure,  synthetic  diamond  manufacturing  methods  must  provide  the  favorable  ther¬ 
modynamic  environment  or  chemical  kinetic  conditions  that  promote  diamond  growth.  The 
three  main  manufacturing  methods  to  create  synthetic  diamond  are  high  pressure,  high  tem¬ 
perature  synthesis  (HPHT);  chemical  vapor  deposition  (CVD);  and  detonation  techniques  to 
create  nanodiamonds.  The  choice  of  diamond  type  for  NV  Quantum  Information  Processing 
depends  on  the  desired  application.  As  nanodiamonds  with  NV  centers  fall  outside  the  scope 
of  this  thesis,  their  synthesis  and  properties  will  not  be  discussed  here. 

HPHT  synthesis  takes  the  thermodynamic  approach  and  provides  the  temperatures  and 
pressures  necessary  for  stable  diamond  formation.  This  method,  which  mimics  the  condi¬ 
tions  by  which  natural  diamonds  are  formed  in  the  Earth’s  mantle,  was  first  demonstrated 
by  Bundy  et  al.  in  1955  at  General  Electric  [39J.  Seed  crystals  (for  single  crystal  diamond 
growth,  the  seed  is  usually  a  seed  diamond)  are  placed  alongside  a  high-purity  graphitic 
carbon  source  and  metal  catalyst  in  an  anvil  press,  at  pressure  around  5  GPa  with  tempera¬ 
tures  exceeding  1400°C.  The  metal  catalyst  (usually  a  solvent  of  transition  metals  like  nickel, 
cobalt,  and/or  iron)  dissolves  the  graphitic  carbon  and  transports  it  to  the  seed  crystal  in 
addition  to  lowering  the  temperature  needed  for  diamond  formation  [40j.  A  temperature 
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gradient  is  applied  to  precipitate  carbon  onto  the  seed  crystal,  which  then  grows  into  the 
desired  diamond. 

CVD,  in  contrast  to  HPHT  synthesis,  creates  chemically  kinetic  conditions  that  favor  the 
growth  of  diamond  over  graphite.  This  allows  lower  temperatures  and  pressures  to  be  used 
during  growth.  A  substrate  (again,  usually  a  diamond  for  single  crystal  growth)  is  placed 
in  a  reactor  along  with  gas  containing  carbon  and  hydrogen.  After  the  gas  is  heated  into 
a  plasma  the  ensuing  carbon  radicals  deposit  onto  the  substrate,  which  in  this  metastable 
environment  can  be  in  the  form  of  graphite  (sp2-hybridized)  or  diamond  (sp3-hybridized). 
Early  methods  in  the  1960’s  used  filaments  to  heat  and  dissociate  the  gas,  but  (to  avoid 
introducing  impurities  into  the  diamond)  microwaves  are  primarily  used  today  [41j.  The 
hydrogen  radicals  etch  off  carbon  that  hybridizes  in  the  sp2  form,  allowing  diamond  to  grow 
layer-by-layer  instead  to  form  plates  a  few  millimeters  in  size  S3 

A  consideration  in  both  of  these  methods  is  the  isotopic  content  of  the  carbon  used  to 
form  the  diamond  lattice.  Naturally  occurring  carbon  is  composed  1.1%  of  the  isotope  13C, 
which  has  a  nuclear  spin  /  =  1/2.  The  dipolar  coupling  between  this  spin  and  the  NV 
electronic  spin  is  a  source  of  decoherence  that  limits  the  coherence  times  T2  and  T2*  °f  the 
NV  center.  If  the  isotopic  content  is  controlled  (e.g.  using  isotopically  pure  12C  methane 
gas  in  CVD  synthesis),  the  abundance  of  this  isotope  decreases  and  the  coherence  time  is 
extended  [38,  43j.  At  the  time  of  this  thesis,  only  naturally  abundant  diamonds  (1.1%  13 C) 
were  available  for  our  study. 

Although  HPHT  synthesis  enables  mass-production  of  diamond  that  can  be  used  for 
numerous  industrial  applications,  for  NV  magnetometry  diamonds  created  by  CVD  are  pre¬ 
ferred  for  the  following  reasons.  First,  as  discussed  below,  a  high  concentration  of  para¬ 
magnetic  substitutional  nitrogen  lowers  the  coherence  time  of  the  NV  center  [44] .  HPHT 
synthesis  typically  produces  diamonds  with  100  —  200  ppm  of  nitrogen  (which,  since  nitrogen 
absorbs  mostly  blue  and  green  light,  gives  the  diamond  a  yellow  color),  while  CVD  produces 
a  range  usually  between  1  ppb  to  a  few  or  several  ppm,  where  the  effects  of  other  impurities 
like  13C  can  become  the  dominant  source  of  decoherence.  Second,  CVD  allows  more  pre¬ 
cise  control  over  the  defects  introduced  into  the  diamond  by  means  of  the  gases  used.  This 
allows  higher-purity  diamonds  to  be  produced  and  altered  to  fit  a  particular  technological 
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application. 


2.2.2  Nitrogen  Incorporation 

In  addition  to  the  way  the  host  diamond  is  grown,  the  nitrogen  content  is  a  crucial  factor 
in  creating  NV  centers  in  the  diamond  lattice.  In  both  HPHT  and  CVD  diamonds,  nitro¬ 
gen  is  a  major  impurity  incorporated  during  growth.  In  addition  to  deliberate  addition  of 
nitrogen  (usually  N2)  during  diamond  growth,  the  residual  air  in  the  growth  chamber,  the 
carbon  gas  sources  (usually  methane),  and  the  metal  catalysts  (nickel,  iron,  and  cobalt  are 
commonly  used)  are  just  a  few  sources  that,  if  left  unchecked,  can  cause  the  creation  of 
unwanted  nitrogen  defects  in  the  lattice.  In  both  HPHT  [45j  and  CVD  [46,  47J  diamonds 
single  substitutional  nitrogen  N5  (also  called  the  PI  center)  is  the  main  form  of  the  nitro¬ 
gen  defect.  If  the  nitrogen  is  not  controlled,  HPHT  growth  may  give  rise  to  hundreds  of 
ppm,  which  becomes  the  dominant  limiting  factor  of  the  coherence  time  of  the  NV  center 
[3  0H  08].  Controlled  nitrogen  incorporation  can  be  achieved  by  the  use  of  nitrogen-getters 
like  aluminium  or  boron  [49J  and  injecting  a  controlled  amount  of  nitrogen  gas  into  the 
growth  chamber.  For  CVD  diamond,  the  amount  of  nitrogen-containing  gas  can  be  varied 
to  get  a  desired  nitrogen  concentration,  but  with  great  care  nitrogen  concentrations  lower 
than  100  ppb  can  be  obtained  ESI 

Another  method  to  incorporate  nitrogen  into  the  lattice  is  through  nitrogen  ion  implan¬ 
tation,  usually  with  the  ions  N+  or  Nj .  The  advantages  of  this  method  are  that  it  allows 
a  greater  control  over  the  nitrogen  dosage  and  introduces  the  extra  vacancies  needed  in  the 
diamond  to  form  NV  centers  (after  subsequent  annealing,  as  explained  below).  The  position 
and  depth  can  also  be  controlled  with  great  precision,  making  this  an  advantageous  ap¬ 
proach  for  single  NV  experiments  or  experiments  with  NV  layers.  For  nitrogen  implantation 
the  use  of  15N  can  allow  implanted  nitrogen  to  be  distinguished  from  existing  nitrogen  in 
the  diamond.  The  disadvantages,  however,  make  this  method  less  appealing  for  ensemble 
magnetometry:  the  optical  and  spin  coherence  properties  of  NV  centers  formed  with  these 
nitrogen  ions  is  usually  worse  compared  to  as-grown  NV  centers.  Two  relevant  reasons  are 
the  higher  preference  to  be  in  the  NV°  charge  state  |£QE21  and  the  shorter  coherence  times 
due  to  radiation  damage  caused  by  the  nitrogen  to  the  diamond  lattice  [38,  [53].  For  this 
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reason,  we  used  the  as-grown  nitrogen  in  our  CVD  diamonds  to  make  NV  centers  and  did 
not  pursue  ion  implantation. 

2.2.3  Vacancy  Production 

Regardless  of  how  the  nitrogen  is  added  to  the  diamond,  to  form  nitrogen- vacancy  centers 
the  corresponding  vacancies  need  to  be  present  in  the  lattice.  For  NV  centers  it  is  believed 
m  that  two  nitrogen  atoms  are  required  to  form  a  NV  center:  one  to  pair  with  a  vacancy 
and  another  to  donate  an  electron  to  make  the  NV  center  negatively  charged.  This  sets  the 
ideal  nitrogen  to  NV  conversion  efficiency  to  50%;  were  this  to  be  exceeded  the  native  NV° 
concentration  would  increase  as  the  available  substitutional  nitrogen  Ng  charge  donors  would 
decrease.  Although  for  CVD  diamonds  NV  centers  are  present  after  growth,  the  nitrogen  to 
NV  conversion  efficiency  is  usually  less  than  0.5%  [53]  •  Thus,  to  create  more  vacancies  and 
enhance  the  NV  concentration,  the  diamond  is  usually  irradiated  with  high-energy  particles 
to  displace  carbon  atoms  from  their  lattice  sites.  The  caveat  to  this  method  is  that  additional 
defects  may  be  introduced  into  the  diamond,  giving  rise  to  additional  paramagnetic  defects 
[38,  071  E3  56J  that  can  lower  the  coherence  time  of  the  NV  center. 

To  minimize  unwanted  damage  to  the  diamond  lattice  while  still  increasing  the  nitrogen 
vacancy  concentration,  the  type  and  energy  of  irradiating  particles  can  be  controlled.  As 
discussed  above  nitrogen  ions  are  an  option  to  create  vacancies,  but  other  ions  and  heavy 
particles  can  be  used  such  as  protons  and  neutrons  [36j.  Gamma  rays  are  also  an  option  |57], 
For  NV  magnetometry,  however,  these  high  energy  and  momentum  particles  are  undesirable 
due  to  the  large  amounts  of  damage  they  create.  Irradiation  of  this  type  can  produce 
knock-on  atoms  [58],  which  occurs  when  the  carbon  atom  is  displaced  from  the  lattice  site 
with  enough  energy  to  itself  displace  further  carbon  atoms.  Additionally,  these  heavier 
particles  can  create  localized  lattice  damage  that  cannot  be  annealed  out,  as  described  in 
the  following  section  [59J,  as  well  as  paramagnetic  vacancy  chains  [38J.  Electrons  on  the 
other  hand  are  light  enough  to  avoid  this  type  of  damage  and  are  more  suited  to  creating 
evenly  distributed  monovacancies  [58].  This  is  the  preferred  irradiation  particle  to  create 
NV  ensembles  optimized  for  NV  magnetometry. 

The  energy  of  the  irradiation  electrons  determines  the  depth  of  penetration.  Experimen- 
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tally  vacancy  creation  has  been  observed  for  incident  electrons  with  energies  greater  than 
180  keV  [60j,  which  gives  the  electron  enough  energy  to  collide  with  the  carbon  atoms  and 
transfer  about  35  eV  of  activation  energy.  This  value  closely  matches  theoretical  predictions 
of  [57J  where  monovacancy  production  is  calculated  to  occur  with  electron  energies  of  at  least 
165  keV.  To  avoid  excess  damage  to  the  diamond  lattice,  it  is  desirable  to  keep  the  electron 
energy  low  enough  to  prevent  the  creation  of  vacancy  chains  (but  high  enough  to  penetrate 
through  the  diamond). 

2.2.4  Thermal  Annealing 

The  final  step  in  NV  creation  is  the  diffusion  of  vacancies  to  the  nitrogen  atoms  embedded 
in  the  lattice,  which  is  accomplished  by  high  temperature  annealing  of  the  diamond.  In 
addition  to  alleviating  the  crystal  strain  caused  by  the  irradiation,  vacancies  begin  to  diffuse 
in  the  lattice  at  600°C  161]  and  get  trapped  by  the  nitrogen  atoms.  Annealing  at  even  higher 
temperatures  causes  certain  defects  to  ‘anneal’  out,  but  can  also  cause  the  creation  of  more 
stable  defects  [38,  55j. 


2.3  Diamonds  used  in  this  Thesis 


Section  22]  showed  that  there  are  numerous  ways  to  grow  and  to  process  a  diamond,  which 
results  in  variable  properties  NV  centers.  For  our  purposes  we  want  to  choose  processing 
recipes  that  best  tailor  the  diamond  for  NV  ensemble  magnetometry.  For  the  type  of  dia¬ 
mond  to  use,  CVD  is  likely  preferred  over  HPHT  due  to  the  lower  nitrogen  content,  higher 
control  of  defect  incorporation  and  isotopic  composition,  and  easier  ability  to  control  the 
geometry  of  the  sample  (e.g.  create  a  layer  with  nitrogen  defects  atop  a  purer  diamond 
substrate).  For  this  study  we  chose  commercially  available  CVD  diamonds  with  an  initial 
nitrogen  concentration  of  0.15  ppm  [62 J  and  a  natural  abundance  of  13 C  (1.1%).  Although 
a  higher  concentration  of  nitrogen  and  a  high  purity  12C  diamond  would  be  better  suited  to 
NV  ensemble  magnetometry,  at  the  time  of  this  study  these  were  the  commercially  available 
CVD  diamonds,  in  addition  to  HPHT  diamonds  at  a  few  hundred  ppm  of  nitrogen.  How¬ 
ever,  the  commercial  availability  gave  us  the  means  to  acquire  a  large  set  of  diamonds  for  a 
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systematic  study  of  diamond  processing,  with  the  intent  of  studying  the  effects  of  irradiation 
and  annealing  and  extracting  results  that  can  be  applied  to  diamonds  suited  for  magnetome- 
try  when  they  become  more  available.  Electron  irradiation  at  a  commercial  company  (Prism 
Gem)  using  1  ME  electrons  is  used  to  produce  vacancies  in  the  diamonds. 


The  diamonds  in  the  study  are  commercially  available  standard  grade  diamond  available 
from  Element  Six.  The  size  of  the  diamonds  are  thin  square  wafers  of  side-length  2.6  mm 
and  thickness  of  0.3  mm  with  un-polished  sides  and  (100)  polished  surfaces.  The  nitrogen- 
14  concentration  listed  by  Element  Six  is  0.15  ppm  (w  2.64  x  1018  cm-3).  The  natural 
isotopic  abundance  of  carbon  is  98.9%  12C  and  1.1%  13C.  Prism  Gem  in  New  York  irradiated 
and  annealed  the  diamonds.  To  ensure  the  radiation  damage  was  uniform  throughout  the 
diamond,  the  1  ME  electron  beam  used  was  1  cm  in  diameter.  A  given  electron  irradiation 
dose  in  electrons  per  cm2  depended  on  the  duration  of  the  irradiation  treatment.  The 
samples  were  placed  on  water  cooled  plates  not  exceeding  66°C,  and  Prism  Gem  predicted 
the  diamond  sample  temperature  was  under  120°C  during  irradiation.  The  samples  are 
flipped  halfway  throughout  the  irradiation.  The  annealing  recipes  all  involved  first  heating 
the  diamonds  to  850°C  in  a  reduced  nitrogen  atmosphere  of  1.0  atm.  A  thermocouple  kept 
track  of  the  temperature  as  the  diamonds  annealed  for  an  hour.  Those  diamonds  slated 
to  receive  a  1250°C  anneal  were  then  moved  to  another  furnace  in  a  100  torr  (0.13  atm) 
hydrogen  atmosphere  for  an  hour.  The  temperature  in  this  chamber  was  recorded  by  a 
pyrometer  pilot  chip.  |Figure  2%2  shows  the  effects  of  annealing  on  diamond  color  for  three 
of  our  samples  that  had  the  same  irradiaton  dose. 


All  twenty-eight  diamonds  we  characterized  are  compiled  in  |Table  2.1j  With  the  excep¬ 
tion  of  three  diamonds  (two  of  them  being  control  diamonds  that  received  no  post-growth 
processing),  all  were  subjected  to  an  initial  anneal  at  850°C  for  one  hour.  After  a  given  dose 
of  irradiation  the  diamonds  were  either  not  annealed,  annealed  at  850°C  for  one  hour,  or 
annealed  850°C  for  one  hour  followed  by  an  anneal  at  1250°C  for  one  hour.  The  range  of 
electron  doses  was  chosen  such  that  the  range  of  vacancy  concentration  before  annealing  be 
less  than,  comparable  to,  and  greater  than  the  nitrogen  concentration  in  the  diamond. 
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Sample 

Pre-anneal  (°C,  1  hr) 

Irradiation  (1 

ME  e“/cm2) 

Post-anneal  (°C,  1  hr) 

Meitner 

none 

none 

none 

Heaviside 

none 

none 

none 

Hyde 

850 

none 

850 

Ayrton 

850 

1  x  1015 

850 

Curie 

850 

5  x  1015 

850 

Leavitt 

850 

1  x  1016 

850 

Robbins 

850 

1  x  1016 

850 

Einstein 

850 

5  x  1016 

850 

Blodgett 

850 

1  x  1017 

850 

Wegener 

850 

5  x  1017 

850 

Bascom 

850 

5  x  1017 

850 

Ramon  y  Cajal 

850 

1.6  x  1018 

850 

Evans 

850 

5  x  1018 

850 

Ladd-Franklin 

none 

5  x  1018 

850 

Thomson 

850 

1.6  x  1019 

850 

Cannon 

850 

5  x  1019 

850 

Diesel 

850 

1  x  1015 

1250 

Hilbert 

850 

5  x  1015 

1250 

Rutherford 

850 

1  x  1016 

1250 

Blau 

850 

X 

o 

05 

1250 

Noddack 

850 

5  x  1016 

1250 

Cori 

850 

1  x  1017 

1250 

Bohr 

850 

5  x  1017 

1250 

Mendeleev 

850 

5  x  1017 

1250 

Stevens 

850 

1.6  x  1018 

1250 

Lorentz 

850 

5  x  1018 

1250 

Anderson 

850 

1.6  x  1019 

1250 

Planck 

850 

5  x  1019 

1250 

Table  2.1:  Diamonds  processing  conditions  of  irradiation  dose  and  annealing  temperature. 
For  sanity,  bookkeeping,  and  a  little  amusement,  diamonds  were  named  after  scientists. 
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1  MeV  Electron-Irradiated 
5x1 018  e/cm2 
CVD  Diamond 


No  850  °C  1250  °C 

Anneal  Anneal  Anneal 


Figure  2-2:  Three  Element  Six  standard  grade  diamonds  used  in  this  thesis  with  the  same 
electron  dose  subjected  to  varying  annealing  temperatures.  The  change  in  color  is  due  to 
the  annealing  out  of  different  defects  at  higher  temperatures. 

2.4  Conclusion 


In  this  chapter  and  the  previous  one  we  have  set  the  stage  for  the  experiments  performed  in 
this  thesis.  IChapterT]  discussed  the  sensitivity  equation,  a  key  metric  in  NV  magnetometery, 
in  terms  of  its  main  variables  and  how  each  can  be  optimized.  In  this  chapter  we  then  looked 
at  how  diamond  processing  techniques  could  be  used  to  create  and  tailor  a  NV  ensemble  in 
diamond  for  the  purpose  of  high-sensitivity  magnetometry. 

The  remainder  of  this  thesis  is  dedicated  to  two  sets  of  measurements  designed  to  test  each 
diamond  for  properties  relevant  to  magnetometry.  [Chapter^] will  describe  photoluminescence 
spectroscopy  measurements  used  to  examine  the  charge  state  fluorescence  ratio  of  NV-  to 


NVU  along  with  the  total  fluorescence  of  the  NV  ensemble.  Chapter  4  details  the  coherence¬ 


time  measurements  that  extract  the  spin  decoherence  time  T2  for  each  diamond.  Finally,  in 


Chapter  5]  the  results  will  be  analyzed  to  obtain  an  optimal  parameter  space  for  diamond 
processing  with  commercially  available  CVD  diamonds. 


50 


Chapter  3 


Charge  State  Fluorescence 
Characterization 

3.1  Overview 

As  outlined  in  the  previous  two  chapters,  the  goal  of  this  thesis  is  to  understand  the  effects  of 
irradiation  and  annealing  on  NV-  formation  and  the  spin  properties  in  diamond  for  ensemble 
magnetometry.  To  accomplish  this  goal  we  fix  the  type  and  energy  of  the  irradiation  particle 
(1  ME  electrons)  and  initial  nitrogen  concentration  (0.15  ppm)  and  systematically  vary  the 
irradiation  dose  and  annealing  temperature.  We  determine  which  combination  optimizes 
the  sample  set  for  magnetometry.  To  quantify  magnetometry  performance,  we  look  at  the 
sensitivity  equation  (77  ~  l/a^/fiT),  which  depends  on  the  fluorescence  contrast  a,  photons 
collected  (3  per  measurement,  and  interrogation  time  T  [7J  30,  32 j  (which  is  bounded  by  the 
NV  coherence  time  T2).  In  this  chapter  we  start  the  characterization  of  the  28  diamonds  used 
in  this  study  by  evaluating  the  contrast  a  and  the  photons  collected  (3  by  photoluminescence 
spectroscopy  measurements.  The  other  variable  in  the  ever-present  sensitivity  equation,  T, 
will  be  addressed  in  the  next  chapter. 

After  discussing  the  role  of  the  contrast  a  and  collected  photons  /3  in  magnetometry 
measurements,  this  chapter  will  describe  the  acquisition  of  photoluminescence  (PL)  spectra 
and  subsequent  analysis  by  fitting  to  a  linear  combination  of  NV°  and  NV-  spectra.  Two 
relevant  quantities  are  extracted  from  the  PL  spectra:  the  total  NV  fluorescence  (which 
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measures  the  collected  photons  f3)  and  the  charge  state  fluorescence  ratio  (which  affects  a). 
The  behavior  of  these  two  quantities  with  irradiation  dose  and  laser  excitation  power  will  be 
shown. 


3.2  Charge  State  Dynamics  and  Magnetometry 


InjChapter  l|we  defined  the  contrast  a  in  terms  of  the  two  signals  acquired  in  a  magnetometry 
measurement,  from  which  the  NV-  population  difference  (and  thus  the  magnetic  field)  is 
calculated.  For  the  two  fluorescence  rate  signals  mo  and  mi,  the  contrast  is  a  =  (mo  — 
mi) /(mo  +  mi).  The  number  of  photons  collected  in  a  measurement  was  similarly  defined  in 
terms  of  mo  and  mp  /?  =  (mo  +  mi) /2tgate,  where  tgate  is  the  fluorescence  collection  time.  To 
optimize  a  and  /3  for  the  magnetic  sensitivity,  a  high  conversion  efficiency  of  a  fixed  amount 


of  single  substitutional  nitrogen  Ns  to  NV-  centers  is  required.  As  mentioned  in  [Chapter 
[2j  to  form  an  NV-  center,  typically  two  nitrogen  atoms  in  the  diamond  lattice  are  required, 
one  to  capture  the  vacancy  and  another  nitrogen  (or  charge  donor)  to  donate  an  electron 
to  give  the  NV  a  negative  charge  [63j.  As  the  initial  Ns  to  NV-  conversion  efficiency  in 
as-grown  synthetic  diamond  is  low  [54J,  irradiation  of  high  energy  particles  or  ions  is  used 
to  create  more  vacancies  that,  upon  subsequent  annealing,  increase  the  NV-  concentration. 
This  process  increases  the  number  of  collected  photons  in  the  signals  mo  and  mi  and  thus 
improves  (3. 

An  unwanted  consequence  of  this  processing  procedure,  however,  is  the  increased  presence 
of  the  neutrally-charged  nitrogen- vacancy  center  NV°  m,  which  does  not  possess  the  desired 
optical  and  spin  properties  needed  for  magnetometry.  If  a  vacancy  captured  by  a  nitrogen 
atom  is  not  located  close  to  an  Ng  atom,  it  will  not  receive  the  needed  electron  to  become 
NV-.  If  the  concentration  of  created  vacancies  from  irradiation  approaches  the  concentration 
of  nitrogen,  NV°  wifi  be  the  dominant  NV  charge  state  in  the  diamond  lattice  as  there 


will  be  fewer  Ns  to  donate  electrons  [23j.  Figure  3-1  shows  the  PL  spectrum  for  Diesel 
(1  x  1015  e-/cm2,  1250°C  anneal)  excited  with  10  mW  of  laser  power.  The  PL  spectrum 
shows  the  zero-phonon-lines  (ZPL)  of  both  the  NV°  (575  nm)  and  NV-  (637  nm)  charge 
states,  as  well  as  their  respective  phonon-sidebands.  The  Raman  line  at  572  nm,  which  arises 
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2000 


Figure  3-1:  Example  photoluminescence  spectra  for  Diesel  (1  x  1015  e-/cm2,  1250°C  anneal) 
at  10  mW  of  laser  excitation.  The  PL  signatures  from  the  NV°  and  NV-  charge  states  are 
also  shown  to  highlight  the  contribution  of  each  charge  state,  as  well  as  the  overlap  of  the 
two  spectra. 


from  inelastic  light  scattering  of  the  laser  photons  with  the  diamond  lattice,  is  also  shown. 
The  overlap  of  the  NV°  PL  spectrum  adds  background  fluorescence  to  NV-  state  readout. 
As  NV°  does  not  respond  to  the  MW  excitation  that  creates  a  population  difference  in  NV-, 
the  overall  photons  collected  /3  increases,  but  the  contrast  a  decreases. 


It  is  thus  desirable  to  determine  irradiation  dose  required  for  a  given  initial  concentration 
of  nitrogen  that  maximizes  the  to  NV-  conversion  efficiency  while  minimizing  the  degrad¬ 
ing  effects  of  NV°,  both  native  NV°  and  that  produced  in  readout.  Even  if  the  conversion 
efficiency  is  high  (low  native  NV°  concentration)  post-processing,  there  are  other  physical 
mechanisms  that  lead  to  the  presence  of  NV°  in  the  lattice  while  under  optical  excitation. 


Recall  from  Chapter  1  that  the  ionization  of  NV  into  NV°  under  optical  excitation  is  a 
two-photon  process,  as  is  the  recombination  process  of  NV°  to  NV-.  Under  the  common 
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operating  conditions  of  low-intensity  green  laser  excitation,  the  dynamics  of  these  competing 
processes  give  rise  to  the  NV  center  on  average  residing  in  the  neutrally-charged  state  30% 
of  the  time  [19,  26j.  In  addition,  experiments  have  shown  that  increasing  the  intensity  of 
optical  excitation  further  photo-ionizes  the  NV  center  to  spend  a  larger  fraction  of  time  in 
the  neutral  charge  state  mm,  since  it  is  a  two-photon  process  as  explained  in 

Given  these  considerations,  we  characterize  the  charge  state  fluorescence  of  our  diamonds 
as  a  function  of  excitation  power  to  provide  a  measure  of  the  collected  photons  /3  and  contrast 
a.  For  the  collected  photons  we  monitor  the  integrated  fluorescence  intensity  of  each  charge 
state  as  a  function  of  irradiation  dose  and  optical  excitation  power  (the  photons  collected  (3 
also  depends  on  the  gate  time  of  fluorescence  collection,  but  this  will  not  be  addressed  in  this 
thesis).  We  provide  a  simple  model  for  explain  the  fluorescence  behavior  with  irradiation 
to  investigate  the  conversion  of  irradiation  dose  to  NV  centers.  Similarly,  for  contrast  we 
measure  the  ratio  of  the  fluorescence  of  each  charge  state  as  a  function  of  irradiation  dose 
and  optical  excitation  power. 


Chapter  1 


3.3  Measurement  Procedure 

To  measure  the  PL  spectra  of  each  diamond,  off-resonant  532  nm  excitation  is  focused  with 
an  objective  to  a  spot  size  of  3.1  /rm  in  diameter  and  depth  of  focus  (Rayleigh  Range)  of 
zr  =  32  /j,m.  A  full  description  of  the  experimental  setup  is  in  |Appendix  A[  The  same 
objective  collects  the  resulting  fluorescence,  which  propagates  to  a  10-/um  multimode  fiber 
that  acts  as  a  pinhole.  This  fiber  directs  the  light  to  a  spectrometer  with  a  CCD  camera  in 
the  readout  plane.  All  frames  taken  by  the  spectrometer/CCD  for  a  single  measurement  are 
averaged  together,  followed  by  normalization  by  the  exposure  time  of  the  measurement.  To 
analyze  the  fluorescence  from  each  charge  state  for  each  diamond  at  different  laser  powers, 
we  fit  each  measured  PL  spectrum  to  a  linear  combination  of  NV°  and  NV-  spectra  [65J,  as 
described  in  |Appendix  B|  The  ratio  of  the  coefficients  forms  the  charge  state  fluorescence 
ratio,  while  the  integral  of  the  coefficient-weighted  basis  spectra  gives  the  fluorescence  of 
each  charge  state. 

Frequently  in  the  literature  the  fluorescence  of  NV-  centers  is  commonly  used  to  estimate 
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their  concentration  in  a  given  diamond,  where  the  integrated  ZPL  fluorescence  from  an 
ensemble  PL  spectra  is  normalized  to  the  integrated  ZPL  of  a  single  NV  center/ensemble 
of  known  concentration  [35,  53,  66j.  The  proper  conversion  of  fluorescence  to  concentration 
is  not  as  trivial  as  it  first  appears,  however.  First,  the  fluorescence  must  be  taken  in  the 
low  laser  intensity  regime,  where  the  fluorescence  is  proportional  to  the  concentration  of 
the  emitters.  At  higher  excitation,  charge  state  conversion  mm  and  saturation  effects 
introduce  a  more  complex  relation  between  the  fluorescence  and  concentration.  Second, 
because  the  NV  center  is  dynamically  switching  between  different  charge  states  under  optical 
excitation,  the  ensemble  fluorescence  measured  is  a  time-averaged  value  [19j.  This  means 
the  charge  state  fluorescence  ratio  is  not  equal  to  the  charge  state  concentration  ratio.  A 
full  rate-equation  analysis  of  the  different  levels  of  NV°  and  NV-  is  required  for  an  accurate 
comparison;  however,  the  full  electronic  structure  (and  associated  transition  rates)  of  NV° 
is  still  the  subject  of  active  investigation.  For  these  reasons  we  report  the  charge  state 
fluorescence  only.  This  is  also  the  most  important  parameter  because  in  practice  this  ratio 
determines  the  contrast  a  and  collected  photons  f3  in  a  magnetometry  measurement. 

We  now  systematically  show  the  results  of  charge  state  fluorescence  measurements,  both 
as  a  function  of  irradiation  dose  and  optical  excitation  power. 


3.4  Charge  State  Fluorescence  Measurements 

3.4.1  Total  Fluorescence  vs.  Irradiation  Dose 


We  first  characterize  the  collected  photons  f3  for  each  diamond  by  calculating  the  total  NV 
fluorescence.  The  total  NV  fluorescence  is  obtained  by  adding  the  fluorescence  from  each 


charge  state,  determined  from  the  linear  combination  of  the  NV  basis  functions.  Figure  3-2 


shows  the  total  NV  fluorescence  vs.  irradiation  dose  D  over  four  orders  of  magnitude  of 
laser  power.  For  clarity  the  data  from  two  annealing  temperatures  are  plotted  separately. 
For  both  annealing  temperatures  and  all  laser  excitation  powers,  we  observe  three  regions  of 
different  behavior:  at  low  electron  irradiation  doses  (D  <  lx  1018  e“/cm2),  the  fluorescence 
steadily  increases;  at  intermediate  doses  (1  x  1018  e“/cm2  <  D  <  1  x  1019  e“/cm2)  the 
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850  °C  Anneal 


1250°C  Anneal 


Figure  3-2:  Total  NV  fluorescence  vs.  irradiation  dose  over  four  orders  of  magnitude  for  dia¬ 
monds  that  underwent  (left)  850°C  annealing  treatment  (right)  1250°C  annealing  treatment. 


fluorescence  saturates,  and  at  the  highest  doses  (D  >  1  x  1019  e-/cm2),  the  fluorescence 
begins  to  decrease.  The  first  two  dose  regions  are  consistent  with  the  discussion  of  vacancy 
production  in  the  lattice  in  |Section  T2|  Higher  irradiation  doses  introduce  more  vacancies 
into  the  lattice,  which  after  annealing  lead  to  the  formation  of  more  NV  centers  (NV- 
and  NV°).  As  the  vacancy  concentration  surpasses  that  of  nitrogen  in  the  diamond,  the 
fluorescence  saturates  as  the  production  of  NV  ceases.  The  behavior  at  the  highest  doses 
where  the  fluorescence  begins  to  decrease  is  not  expected  from  considering  nitrogen  and 
vacancy  concentrations  alone.  One  possible  cause  could  be  the  creation  of  vacancy-related 
defects  that  quench  the  concentration  of  NV  centers  [55j. 

To  compare  the  affects  of  annealing  temperature,  we  look  at  both  temperatures  simulta¬ 
neously.  |Figure  3^31  displays  the  total  fluorescence  of  all  diamonds  as  a  function  of  irradiation 
dose  for  low  (100  /iW)  and  high  (10  mW)  laser  excitation.  In  terms  of  the  behavior  discussed 
above,  the  fluorescence  of  1250°C  annealed  diamonds  saturate  at  a  lower  dose  (~  1  x  1017 
e-/cm2)  than  that  of  the  850°C  annealed  diamonds  (~  5  x  1017  e-/cm2).  This  means  the 
efficiency  of  NV  formation  for  a  given  dose  is  higher  for  1250°C  anneal  diamonds.  This  can 
be  understood  in  terms  of  vacancy  diffusion:  in  a  diamond  subjected  to  a  hotter  annealing 
temperature,  the  vacancies  diffuse  farther  in  that  lattice  for  a  given  amount  of  time,  giving 
them  a  higher  probability  to  be  trapped  by  a  nitrogen  atom.  Note,  however,  that  at  all 
powers  shown  the  highest  fluorescence  occurs  at  850°C,  the  decrease  of  fluorescence  at  the 
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highest  doses  seems  to  seems  to  be  independent  of  the  annealing  temperature,  however,  since 
it  occurs  at  about  the  same  dose  (~  1  x  1019  e"/cm2)  for  both  temperatures. 


Irradiation  Dose  (e'/cm2) 


Irradiation  Dose  (e'/cm2) 


Figure  3-3:  Total  NV  fluorescence  vs.  irradiation  dose  for  all  diamonds  at  (left)  100  /iW  of 
laser  excitation  before  the  objective  and  (right)  10  mW. 


3.4.2  Fluorescence  Model 

To  further  investigate  the  effect  of  annealing  temperature  on  the  conversion  efficiency  of 
irradiation  electrons  to  NV  centers,  we  fit  the  above  fluorescence  data  to  a  simple  model. 
We  ignore  the  influence  of  other  charge  traps  in  the  diamond,  such  as  lattice  dislocations,  and 
focus  on  the  vacancy  concentration  created  by  a  given  irradiation  dose.  Electron  irradiation 
creates  vacancies  V°  in  the  diamond  which,  when  trapped  by  a  single  substitutional  nitrogen 
N^,  form  NV°  [63j : 

Ng  +  V°  ->•  NV°. 

Interaction  with  a  nearby  induces  charge  transfer  to  form  NV-: 

Ng  +  NV°  ->•  N+  +  NV-. 

We  assume  three  dose  ranges  in  this  model:  at  low  irradiation  dose,  the  additional  NV  centers 
(both  charge  states)  formed  are  negligible  compared  to  the  initial  NV  concentration  present 
after  diamond  synthesis.  At  intermediate  doses,  when  the  vacancy  concentration  created  is 
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comparable  to  the  nitrogen  concentration,  the  growth  in  NV  concentration  is  approximately 
linear.  At  higher  doses,  when  the  vacancies  created  far  exceed  the  nitrogen  content  in  the 
diamond,  the  NV  concentration  saturates.  We  ignore  the  fluorescence  at  doses  higher  than 
(5  x  1018  e-/cm2),  as  it  falls  outside  the  scope  of  this  model  and  is  not  needed  to  investigate 
irradiation  to  NV  conversion.  We  also  assume  that  the  concentration  of  NV  centers  is  linearly 
proportional  to  the  fluorescence,  which  is  valid  at  laser  excitation  below  saturation. 

The  functional  form  of  the  model  is 


fNv(D)  =  A  ([NVi]  +  [Nt]  (l  -  exp 


D 


(1.76  x  1017cm_3/ppm)[Nr] 


(3.1) 


where  Jnv{D)  is  the  total  NV  fluorescence  as  a  function  of  irradiation  dose  D,  [NV*]  is 
the  initial  concentration  of  NV  centers  in  ppm  pre-irradiation,  [N^]  is  the  total  nitrogen 
concentration  in  ppm  in  the  diamond  pre-irradiation,  A  is  a  scaling  factor  to  convert  from 
concentration  to  fluorescence  (ppm-1),  and  B  is  a  parameter  relating  the  conversion  of 
irradiation  electrons  to  NV  centers  (V°/e-/cm).  To  fit  the  model  we  fix  the  initial  nitrogen 
concentration  at  0.15  ppm,  which  is  consistent  with  the  typical  nitrogen  concentrations  found 
in  commercially  available  CVD  diamond  from  Element  Six  |62j .  Figure  3-4  shows  an  example 


of  the  model  applied  to  diamonds  incident  with  70  mW  optical  excitation.  The  two  highest 
doses,  also  shown  in  the  figure,  were  not  included  in  the  fitting. 


Figure  3-5  depicts  the  fit  parameters  as  a  function  of  excitation  power.  The  initial 


NV  concentration  is  consistently  in  the  range  of  1-10  ppb,  corresponding  to  an  N  to  NV 
conversion  efficiency  between  0.6  -  6%.  The  parameter  A  increases  and  shows  the  beginning 
of  saturation  at  high  excitation  powers,  which  is  expected  since  the  assumption  that  the  NV 
concentration  is  proportional  to  the  fluorescence  is  not  valid  past  saturation  (calculated  to 


be  5  mW  in  this  experimental  setup,  see  Appendix  B).  The  irradiation  to  vacancy  conversion 


factor  B  is  roughly  constant  until  powers  past  4  mW  for  both  annealing  temperatures,  at 
which  they  both  begin  to  decline.  Since  this  parameter  describes  the  conversion  of  dose  to 
vacancy  concentration,  it  should  not  depend  on  laser  power.  The  decline  is  probably  due  the 
breakdown  of  the  assumption  that  the  fluorescence  is  proportional  to  concentration  at  high 
powers. 
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Figure  3-4:  Sample  model  fitting  for  fluorescence  measurements  taken  at  70  mW  excitation 
power.  Red  corresponds  to  850°C  annealing  and  blue  to  1250°C.  Note  that  the  two  unpro¬ 
cessed  diamonds  are  included  in  both  data  sets.  The  diamonds  at  doses  exceeding  5  x  1018 
e“/cm2  are  also  shown,  distinguished  from  the  main  data  set  as  these  were  not  used  for  the 
fitting. 


To  analyze  the  conversion  efficiency  of  irradiation  electron  dose  to  NV  production,  we 
average  the  B  values  for  powers  less  than  6  mW.  For  850°C  we  have  B$5o  =  0.096  ±  0.028 
V°/e-/cm  and  for  1250°C  we  have  -£>1250  =  1-347  ±  0.601  V°/e_/cm.  That  is  about  0.1 
vacancies  are  formed  per  1  ME  electron  per  centimeter  at  850°C,  while  about  1.3  vacancies 
are  created  per  electron  at  1250°.  As  mentioned  above,  one  explanation  for  the  higher  B 
value  is  the  greater  diffusion  vacancies  experience  at  higher  temperatures. 


3.4.3  Total  Fluorescence  vs.  Optical  Power 


The  total  fluorescence  of  each  NV  charge  state  is  plotted  against  optical  excitation  power  for 


850°C  annealing  in  Figure  3-6  and  for  1250°C  annealing  in  figure  Figure  3-7 L  Both  charge 


1Note  that  Bascom  (5  x  1017  e“/cm1 2,  850°C  anneal),  whose  PL  spectra  serves  as  the  NV°  basis  function, 

shows  a  decrease  in  fluorescence  with  higher  power;  this  is  a  systematic  error  that  can  be  fixed  with  a 
refinement  of  the  basis  functions.  For  the  entire  thesis,  data  for  Bascom  at  high  laser  power  deviates  from 
expected  trends  due  to  this  systematic  error. 
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Figure  3-5:  Fitting  parameters  for  the  fluorescence  model  as  a  function  of  laser  power.  The 
initial  nitrogen  concentration  Nt  is  fixed  at  0.15  ppm.  For  all  plots,  the  error  bars  are 
extracted  from  the  68%  confidence  intervals  obtained  from  the  fitting;  those  associated  with 
A  are  encompassed  within  the  points. 


states  show  the  beginning  of  saturation  in  fluorescence  as  the  optical  power  is  increased.  To 
determine  what  the  saturation  power  is,  we  fit  a  simple  model[l9j  the  fluorescence  /  of  each 
charge  state  to  the  form 

m  =  Airrr-  (3-2> 

r  -T  rsat 

where  P  is  the  laser  power  measured  before  the  objective,  Psat  is  the  saturation  power,  and 
A  is  a  fitting  parameter. 

shows  two  examples  of  this  fit.  For  the  diamond  Bohr  (5  x  1017  e"/cm2, 
1250°C  anneal)  the  NV°  saturation  power  is  extracted  to  be  Psat,o  =  160  mW  and  for 
NV-  is  Psat-  =  9.2  mW.  Similarly,  for  the  diamond  Lorentz  (5  x  1018  e“/cm2,  1250°C 
anneal)  Psat,o  =  158  mW  and  Psat-  =  55  mW.  A  significant  number  of  diamonds  do  not 
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Figure  3-6:  NV  fluorescence  vs.  optical  excitation  power  of  all  diamonds  annealed  at  850°C 
for  the  (left)  neutral  charge  state  and  (right)  negative  charge  state. 
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Figure  3-7:  NV  fluorescence  vs.  optical  excitation  power  of  all  diamonds  annealed  at  1250°C 
for  the  (left)  neutral  charge  state  and  (right)  negative  charge  state. 
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fit  the  saturation  model  above,  like  Bohr.  For  diamonds  where  the  saturation  model  was 
successful,  the  saturation  powers  of  NV-  were  all  higher  than  the  calculated  5  mW  for 
this  experimental  setup.  This  suggests  that  saturation  may  not  be  the  only  mechanism 
occurring  at  high  excitation  power.  Further  analysis  into  the  charge  state  dynamics,  and 
perhaps  fluorescence  measurements  on  a  smaller  imaged  fluorescence  region  (see 
0.  is  needed  to  explain  this  discrepancy. 


Appendix 


Figure  3-8:  Saturation  model  fitting  for  (left)  Bohr  (5  x  1017  e-/cm2,  1250°C  anneal)  and 
(right)  Lorentz  (5  x  1018  e“/cm2,  1250°C  anneal).  The  NV-  curve  for  Bohr  is  scaled  for 
comparison  with  NV°. 


3.4.4  Charge  State  Fluorescence  Ratio  vs.  Irradiation  Dose 

We  now  analyze  the  contrast  a  by  looking  at  the  charge  state  fluorescence  ratio  of  NV-  to 

[9]  shows  the  charge  state  fluorescence  ratio  versus  irradiation  dose  D  over  four  orders  of 
magnitude  of  laser  power.  For  clarity  the  data  from  two  annealing  temperatures  are  plotted 
separately.  For  all  laser  powers,  the  ratio  decreases  for  doses  up  to  1  x  1018  e-/cm2.  For 
doses  past  1  x  1018  e“/cm2,  excitation  power  dependence  arises.  At  low  laser  power  the 
ratio  either  continues  to  decrease  or  levels  off;  for  higher  laser  power  the  ratio  increases.  At 
low  doses  the  drop  in  ratio  can  be  attributed  to  the  increase  in  NV°  production  as  more 
vacancies  pair  off  with  Ns,  which  decreases  the  number  of  electron  donors  to  form  NV-. 
The  behavior  at  high  dose  is  not  expected:  similar  to  the  quenching  of  fluorescence  at  these 


NV°.  This  is  obtained  from  the  ratio  of  the  fitting  coefficients  of  NV  to  NV°.  Figure  3- 
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Figure  3-9:  Charge  state  fluorescence  ratio  NV-/NV°  vs.  irradiation  dose  over  four  orders 
of  magnitude  for  diamonds  that  underwent  (left)  850°C  annealing  treatment  (right)  1250°C 
annealing  treatment. 


same  doses  seen  in|Section  3.4. 1[  something  is  causing  the  charge  state  ratio  to  increase.  As 
mentioned  earlier,  the  production  of  a  vacancy-related  defect  could  be  the  cause,  if  it  were 
to  influence  the  charge  switching  dynamics  between  NV-  and  NV°.  This  influence  would 
increase  at  higher  laser  power,  as  seen  in  the  next  section. 


Irradiation  Dose  (e7cm2)  Irradiation  Dose  (eVcm2) 


Figure  3-10:  Charge  state  fluorescence  ratio  NV  /NV°  vs.  irradiation  dose  for  all  diamonds 
at  (left)  100  pW  of  laser  excitation  before  the  objective  and  (right)  10  mW. 


To  see  the  effects  arising  from  the  different  annealing  temperatures,  we  look  at  both 


temperatures  simultaneously.  jFigure  3-10|  shows  the  charge  state  fluorescence  ratio  as  a 
function  of  electron  irradiation  dose  for  all  28  diamonds  at  low  (100  pW)  and  high  (10  mW) 
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excitation  power.  For  doses  below  5  x  1017  e“/cm2,  diamonds  annealed  at  850°C  generally 
had  a  higher  ratio  than  those  annealed  at  1250°C.  This  decrease  can  be  attributed  to  a 
decrease  in  charge  donors  needed  to  create  NV-.  This  agrees  with  the  vacancy  diffusion 
discussion  above,  where  it  was  hypothesized  that  the  vacancies  created  during  irradiation 
diffuse  farther  away  from  their  initial  location,  which  increases  their  probability  to  be  trapped 
by  a  nitrogen  atom.  This  led  to  a  higher  total  NV  fluorescence  at  lower  doses.  For  the  charge 
state  fluorescence  ratio,  the  creation  of  more  NV  centers  at  lower  doses  means  NV°  is  also 
being  created,  which  means  the  ratio  should  be  lower.  Another  explanation  could  be  the 
creation  of  stable  defects  at  this  annealing  temperature  that  are  stable  charge  traps,  robbing 
electrons  needed  for  NV-  to  form. 


3.4.5  Charge  State  Fluorescence  Ratio  vs.  Optical  Power 


We  finish  our  characterization  of  the  charge  state  fluorescence  by  looking  at  the  charge 
state  fluorescence  ratio  as  a  function  of  excitation  power.  The  charge  state  fluorescence 
ratio  is  plotted  as  a  function  of  laser  power  for  both  annealing  temperatures  in  |Figure  3Al 


For  doses  below  1  x  1016  e-/cm2  for  850°C,  5  x  1015  e-/cm2  for  1250°C,  and  for  the  two 
unprocessed  diamonds  (not  shown)  the  ratio  slightly  increases  at  first  but  then  falls  off  with 
higher  excitation  power.  For  doses  up  to  1  x  1018  the  ratio  decreases  monotonically  for  all 
excitation  powers.  This  is  consistent  with  the  literature,  where  increase  laser  intensity  ionizes 
NV-  to  NV°  [19,  El,  64J.  At  doses  above  1  x  1018  e“/cm2,  however,  the  magnitude  of  the 
ratio  begins  to  increase  at  the  higher  laser  power  (but  still  decreases  as  laser  power  increases), 
and  then  for  both  anneal  temperatures  the  two  diamonds  with  doses  above  1  x  1019  e-/cm2 
the  result  is  notable:  the  charge  state  fluorescence  ratio  increases  with  the  laser  excitation. 
Like  the  increase  in  ratio  with  laser  power  seen  in  the  last  section,  this  is  contrary  to  what 


we  expect.  Figure  3-12  shows  the  PL  spectra  of  Cannon  (5  x  1019  e-/cm2,  850°C  anneal) 
at  all  laser  powers  used.  As  the  excitation  power  increases  the  ZPL  of  NV-  is  increasingly 
distinguishable  from  the  NV°  phonon  sideband. 
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Figure  3-11:  Charge  state  fluorescence  ratio  NV“/NV°  vs.  optical  excitation  power  as 
measured  before  the  dichroic  mirror  in  the  setup,  (top)  850°C  annealing  treatment  (bottom) 
1250°C  annealing  treatment. 
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Figure  3-12:  Photoluminesce  spectra  of  the  diamond  Cannon  (5x  1019  e_/cm2,  850°C  anneal) 
at  all  laser  powers.  Contrary  to  diamonds  irradiated  at  a  lower  dose,  the  NV-  contribution 
to  the  PL  spectra  increases  as  a  function  of  laser  power. 
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3.5  Conclusion 


In  this  chapter,  we  have  described  the  role  that  fluorescence  from  NV°  and  NV-  play  in 
NV  magnetometry.  Photoluminescence  spectra  of  each  diamond  were  measured  for  laser 
intensities  spanning  over  four  orders  of  magnitude  and  fit  to  a  linear  combination  of  individual 
charge  state  spectra.  From  this  fitting  we  calculated  the  total  NV  fluorescence  and  charge 
state  fluorescence  ratio.  For  irradiation  doses  below  1  x  1018  e"/cm2  the  behavior  of  each 
quantity  agreed  with  a  simple  picture  of  adding  more  vacancies  to  a  fixed  concentration 
of  nitrogen.  The  total  fluorescence  from  both  charge  states  increases  with  dose  as  more 
NV  centers  are  created  until  saturating  as  the  vacancy  concentration  is  believed  to  exceed 
that  of  nitrogen.  Similarly,  the  charge  state  fluorescence  ratio  decreases  with  dose  as  the 
available  substitutional  nitrogen  electron  donors  are  used  to  form  more  NV°  centers.  For 
different  annealing  temperatures,  a  simple  model  of  fluorescence  vs.  dose  revealed  a  factor 
of  ~  10  increase  in  irradiation  electron  to  NV  conversion  efficiency  for  1250°C  annealed 
diamonds  over  850°  C,  which  is  attributed  to  the  higher  diffusion  rate  of  vacancies  during 
higher  temperatures.  This  allows  more  vacancies  to  be  trapped  by  nitrogen  atoms  for  a  given 
dose;  it  also  explains  the  lower  charge  state  fluorescence  ratio  for  1250°C  at  low  dose. 

At  higher  irradiation  doses  the  data  behave  contrary  to  what  is  expected  from  our  simple 
model.  The  total  fluorescence  is  quenched  with  the  addition  of  more  vacancies,  and  the 
charge  state  fluorescence  ratio  increases.  Even  more  striking  are  the  diamonds  at  1.6  x  1019 
e“/cm2  and  5  x  1019  e"/cm2,  which  show  a  monotonic  increase  in  charge  state  fluorescence 
ratio  under  increased  laser  excitation,  contrary  to  the  expected  ionizing  role  provided  high 
laser  intensities.  It  is  proposed  that  these  trends  can  be  attributed  to  effects  lying  outside 
the  picture  of  only  nitrogen  and  monovacancies  existing  in  the  lattice,  as  the  NV  fluorescence 
would  not  decrease  past  its  saturated  value,  nor  would  the  charge  state  fluorescence  ratio 
increase  with  the  addition  of  more  vacancies  (NV-  would  ideally  not  exist  at  all  in  the 
lattice  in  the  absence  of  optical  illumination).  High  irradiation  doses  are  known  to  give 
rise  to  vacancy  chains  and  other  vacancy- related  defects  [38,  55j.  These  defects  could  act 
as  charge  donors  to  increase  the  NV-  fluorescence  under  illumination.  In  addition,  the 
defects  could  include  nitrogen  atoms  to  form  defect  complexes,  which  would  decrease  the 
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total  concentration  of  NV  centers  and  thus  the  total  fluorescence.  Further  characterization 
of  defects  present  in  the  highly-irradiated  diamonds  is  required. 

This  data  serves  as  the  basis  to  characterize  the  collected  photons  f3  and  fluorescence 
contrast  a  for  each  diamond.  However,  before  we  can  analyze  the  relative  magnetic  sensitivity 
of  each  diamond,  the  coherence  time  scales  must  also  be  measured. 
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Chapter  4 


Coherence  Time  Characterization 


4.1  Overview 

To  characterize  an  ensemble  of  color  centers  for  magnetometry,  knowledge  of  the  coherence 
times  is  required.  By  analogy  with  nuclear  magnetic  resonance  spectroscopy,  ensembles  are 
characterized  by  two  decay  times,  namely  7\  and  T2  m-  The  relaxation  time  Tf  describes 
population  decay  to  the  ground  state,  while  the  dephasing  time  describes  the  randomizing 
of  phase  between  two  eigenstates.  The  dephasing  time  is  separated  into  two  categories:  the 
inhomogeneous  spin-dephasing  time  T|  characterizes  the  influence  of  inhomogeneities  in  the 
environment,  such  as  temporal  fluctuations  in  the  local  magnetic  held  and  variations  in  the 
local  spin  environment  of  each  NV  center.  The  homogeneous  spin-dephasing  time  T2,  or  spin- 
decoherence  time,  characterizes  the  intrinsic  spin  environment  apart  from  the  influence  of 
external  inhomogeneities  (such  as  applied  MW  and  static  Helds).  Static  dephasing  influences 
are  “focused”  out  in  measuring  the  T2  time,  as  explained  below.  Different  magnetometry 
protocols  depend  on  different  coherence  time  scales. 

In  this  chapter,  we  characterize  the  coherence  times  of  NV  ensembles  by  measuring  the 
three  time  scales  T2*,T2,  and  Tf.  The  role  coherence  time  plays  in  the  sensitivity  equation 
and  the  different  AC  and  DC  magnetometry  regimes  will  be  discussed.  The  emphasis  in  this 
thesis  will  be  on  T2  measurements;  preliminary  T2  and  Ti  measurements  will  also  be  shown. 
For  T2  an  empirical  model  of  the  dependence  of  irradiation  dose  will  be  presented. 
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4.2  Coherence  Time  and  Magnetometry 


As  discussed  in  Chapter  2,  magnetic  sensitivity  (r)  &  1  /oiy/fiT)  depends  on  the  fluorescence 
contrast  a,  total  fluorescence  f3,  and  NV  interrogation  time  T.  High-sensitivity  magnetom¬ 
etry  benefits  from  long  coherences,  which  approach  the  sampling  time  of  the  measurement. 
The  NV  ensemble,  however,  can  only  maintain  phase  coherence  for  a  finite  time  before  the 
surrounding  noise  bath  of  the  environment  causes  ensemble  decoherence.  Past  these  time- 
scales  the  contrast  a  of  the  NV  fluorescence  degrades,  yielding  a  net  loss  in  sensitivity. 
Optimizing  these  time  scales  is  thus  critical  to  an  NV  magnetometer. 


The  first  time  scale, which  is  relevant  for  DC  magnetometry,  is  the  inhomogeneous  spin¬ 
dephasing  time  T2.  The  two  most  common  schemes  to  measure  low  frequency  magnetic 
fields,  CAW-ODMR  and  Ramsey,  both  depend  on  this  time  scale:  for  CAW-ODMR  the 
linewidth  of  the  measured  resonances  scales  as  ~  1/T2,  while  for  Ramsey-based  methods, 
the  NV  ensemble  can  only  maintain  phase  coherence  during  the  free  precession  time  at  time 
scales  equal  to  or  less  than  T2*.  The  magnetic  field  sensitivity  in  both  of  the  above  schemes  are 
thus  optimized  when  the  measurement  time  T  is  of  similar  duration  as  the  inhomogeneous 
spin-dephasing  time  T2  [7,  30j.  The  physical  mechanism  for  spin-dephasing  is  a  result  of 
spatial  and  temporal  field  inhomogeneities  experienced  by  each  NV  center.  For  high-purity 
(low  nitrogen  concentration,  less  than  10  ppm  [48|)  diamond  with  a  natural  abundance  of 
13C  (1.1%),  the  dominant  source  of  dephasing  arises  from  the  hyperfine  dipolar  coupling  of 
the  NV  electronic  spin  with  the  13 C  nuclear  spin.  Dipolar  interactions  with  electron  spins 
from  nitrogen  impurities  (including  other  NV  centers)  also  contribute  to  dephasing. 

In  AC  magnetometry  schemes,  the  relevant  time  scale  shifts  from  T2*  to  T2.  Tailored  pulse 
sequences  can  be  used  to  refocus  static/low  frequency  dephasing  influences  in  the  environ¬ 
ment,  which  extends  the  ability  to  maintain  phase  coherence  to  the  longer  spin-decoherence 
time  T2.  For  NV  ensembles,  T2  can  exceed  T2*  by  a  factor  of  100  [68j  Similar  to  the  DC 
sensing  scheme,  an  AC  measurement  time  of  T  ~  T2  yields  the  optimal  sensitivity.  Like  T2 , 
contributors  to  decoherence  are  dipolar  coupling  to  the  13C  nuclei  and  other  paramagnetic 
impurities,  such  as  the  electron  spins  from  single  substitutional  nitrogen. 


Finally,  the  third  time  scale  to  consider  is  the  spin-lattice  relaxation  time  T\.  This  time 
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scale  characterizes  the  time  for  the  NV  ensemble  to  decay  to  the  thermal  equilibrium  state. 
This  is  the  result  of  interactions  with  lattice  phonons  that  cause  spin-flips.  Ultimately  T2  is 
limited  to  2Ti,  but  it  is  beyond  the  scope  of  this  thesis  to  pursue  methods  of  improving  Tj 
in  our  diamonds. 

The  remainder  of  this  chapter  will  discuss  the  various  pulse  sequences  used  to  measure  the 
three  coherence-times  discussed  above,  with  an  emphasis  on  the  spin-decoherence  time  T2. 
Each  section  first  details  the  measurement  procedure  and  then  the  experimentally-acquired 
results.  The  experimental  measurements  taken  before  each  coherence  time  measurement 
(magnetic  bias  held  alignment  to  an  NV  axis,  determination  of  the  resonances,  etc.)  is  found 
in|Appendix  C[ 


4.3  Spin-Dephasing  Time  T2*  Measurements 


4.3.1  Measurement  Procedure 


The  method  to  measure  T2  is  based  on  the  Ramsey  pulse  sequence  (see  Figure  4-1  (top)[), 
where  after  optical  pumping  to  the  |0)  electronic  ground  state,  |0),  a  microwave  (MW)  7r/2 
pulse  (calibrated  by  a  Rabi  measurement,  see  Appendix  C[)  projects  the  NV  spin  into  a 
superposition  of  |0)  and  |  —  1)  (where  a  static  held  B0  is  applied  to  split  the  |  ±  1)  states). 
The  NV  spin  ensemble  precesses  about  the  static  held  during  a  free  precession  time  r.  A 
second  MW  tt/2  pulse  projects  the  spin  back  into  the  { |0) ,  |  —  1)}  basis  (or  manifold)  where 
the  single-photon  counter  (SPCM)  measures  the  huorescence  m i  with  a  gate  pulse.  The 
sequence  is  repeated  again,  except  a  MW  3tt/2  pulse  is  applied  at  the  end  instead  of  a  MW 
7t/2  pulse  and  huorescence  mo  is  recorded.  The  two  signals  are  subtracted  and  normalized 
to  form  the  signal  s  =  (mi  —  m0)/(mi  +  m0)  to  cancel  out  common-mode  noise. 

The  envelope  of  the  signal  s  is  modeled  as  a  stretched  exponential  with  stretching  pa¬ 
rameter  p: 

3 

s(r)  =  exp[(— t/T2*)p]  ^cos(2tt/;t  -  fc),  (4.1) 


2=1 


where  r  is  the  free  precession  time,  /,  are  the  detunings  from  the  three  NV  hyperhne  res¬ 
onances  that  were  determined  by  a  Fast-Fourier  Transform  of  the  Ramsey  fringes,  <pt  are 


71 


Optical 


Gate 


MW 


Pump  Detection  !  Pump  Detection 


m1  \ 

i 

m0 

T 

I  ;  hh 

n  n  \  7r  : 

In 

2  2  2  2 


Free  Precession  Time  r  (ns) 

Figure  4-1:  (top)  Ramsey  pulse  sequence  used  to  measure  spin-dephasing  time  T2 .  (bottom) 
Example  Ramsey  curve  from  the  diamond  Rutherford  (1  x  1016  e_/cm2,  1250°C  anneal), 
where  the  calculated  spin-dephasing  time  is  T2*  =  719  ±  41.6  ns. 


phase  shifts,  and  p  is  the  stretched  exponential  factor.  Figure  4-1  (bottom) |  shows  a  sample 
measured  signal.  The  detuning  of  the  MW  field  from  the  three  hyperhne  resonances  creates 
the  beating  modulation  of  the  signal.  To  observe  enough  Ramsey  fringes  before  dephasing 
diminishes  the  signal,  the  MW  drive  was  detuned  by  8  MHz  from  the  central  NV  hyperhne 
resonance,  with  the  shortest  detuning  (6  MHz)  having  a  period  of  about  167  ns. 


4.3.2  Results  and  Limitations 

Coherence  measurements  were  carried  out  on  the  matrix  of  irradiated  diamond  in  ITable  2.H 


Figure  4-2  shows  the  measured  T2*  values.  Note  that  Mendeleev  (5  x  1017  e“/cm2,  1250°C 
anneal)  and  Cannon  (5  x  1019  e“/cm2,  850°C  anneal)  were  excluded  due  to  insufficient 
signal-to-noise.  Error  bars  are  extracted  from  the  68%  confidence  intervals  obtained  from 
the  dephasing  curve  fitting.  These  values  lower  bounds  to  the  intrinsic  T2*  of  the  diamond, 
which  for  dilute  nitrogen  concentrations  is  limited  mainly  by  the  dipolar  interactions  of  13C  to 
be  about  1  ps  [T].  Preliminary  measurements  did  not  address  magnetic  held  inhomogeneity 
and  temperature  huctuations,  which  may  limit  the  effective  spin-dephasing  times  reported. 
Further  investigation  is  needed  to  determine  if  the  T2*  values  are  limited  by  either  irradiation 
and  annealing  effects  or  systematic  effects. 

Other  methods  exist  to  measure  T2  .  (1)  Pulsed-ODMR  operates  similar  to  CAW-ODMR 
(see  Appendix  C)  by  measuring  the  huorescence  of  the  NV  ensemble  at  or  near  resonance, 


but  the  laser  is  switched  off  during  pulsed  MW  excitation.  By  reducing  the  MW  power  and 
adjusting  the  MW  ir  pulse,  the  natural  linewidth  Av  can  be  extracted  from  the  huorescence 
lineshape.  Avoiding  laser  power  broadening,  the  spin-dephasing  time  is  then  T2  =  l/(itAv) 
130J.  Pulsed-ODMR  (like  Ramsey)  has  the  advantage  of  eliminating  laser  power  broadening, 
which  allows  the  laser  excitation  to  be  set  high  to  increase  the  huorescence  signal.  This 
method  can  also  be  utilized  through  CAW-ODMR.  (2)  Another  method  to  measure  T2  takes 
advantage  of  the  full  spin-1  nature  of  the  NV  center  and  observe  quantum  beating  between 
bright  and  dark  superposition  states  of  |  ±  1)  magnetic  sublevels  [69J.  The  advantage  of  this 
method  is  to  eliminate  the  role  of  temperature  and  strain-dependent  zero-held  frequency 
[70|  by  measuring  a  beat  frequency  dehned  only  in  terms  of  the  external  magnetic  held 
and  fundamental  constants.  Realizing  this  for  NV  ensembles  is  more  challenging,  however, 
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Irradiation  Dose  (eVcm2) 


Figure  4-2:  Spin-dephasing  times  T|  as  a  function  of  dose  measured  with  the  Ramsey  pulse 
sequence.  Due  to  environmental  factors  that  were  not  explicitly  controlled,  these  values  are 
to  be  seen  as  a  lower  bound  of  the  intrinsic  T2*  of  the  diamond,  not  limited  by  external 
environmental  factors,  the  dotted  line  indicate  the  expected  13C-limited  T2  =  1  /rs  [lj. 

because  only  one  NV  axis  can  be  addressed  at  a  time  and  both  |1)  <->■  |0)  and  |  —  1)  -f-»  |0) 
resonances  need  to  be  driven  at  the  same  time. 


4.4  Spin-Decoherence  Time  T 2  Measurements 


4.4.1  Measurement  Procedure 


The  Spin-echo  (also  known  as  Hahn  Echo)  sequence  (see  Figure  4-3  (top)[),  which  is  used 
to  measure  T2l  builds  on  the  Ramsey  pulse  sequence  by  inserting  an  additional  MW  7r 
pulse  between  the  ir/2  pulses.  As  in  the  Ramsey  sequence,  after  optical  pumping  to  the  |0) 
electronic  spin  state  the  NV  spin  ensemble  is  prepared  into  a  superposition  of  the  |0)  and 
|  —  1)  states  with  a  7t/2  pulse.  The  spin  ensemble  then  acquires  phase  from  the  environment 
during  a  free  precession  time  of  length  r/2.  An  additional  MW  ir  pulse  is  applied  and 
switches  the  sign  of  the  phase  accumulation  for  the  second  free  precession  time.  For  static 
magnetic  fields,  this  means  that  any  phase  accumulated  before  the  MW  ir  pulse  is  “echoed” 
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Figure  4-3:  (top)  Pulse  sequence  used  to  measure  the  spin-decoherence  time  T2-  (bottom) 
Example  signal  from  the  diamond  Noddack  (5  x  1016  e“/cm2,  1250°C  anneal),  where  the 
calculated  spin-decoherence  time  is  T2  =  450  ±2.59  /js .  The  collapses  and  revivals  in  the 
signal  are  due  to  the  Larmor  precession  of  13 C  nuclei  in  the  surrounding  environment. 


out  after.  Phase  accumulation  from  local  magnetic  held  fluctuations  is  also  cancelled  out 
(‘refocused’)  assuming  they  do  not  vary  appreciably  on  the  time  scale  of  the  free  precession. 
Precessing  again  for  a  time  r/2,  a  second  MW  7r/2  pulse  rotates  the  spin  ensemble  back 
into  the  { |0) ,  |  —  1)}  basis  and  a  gate  pulse  from  the  SPCM  records  fluorescence  m\.  As 
in  Ramsey  measurements,  the  readout  pulses  alternates  between  7t/2  and  37t/2  and  output 
signals  from  tuq  and  m\  are  subtracted  to  remove  common-mode  noise. 

Figure  4-3  (bottom) |  shows  a  sample  spin-echo  signal.  The  collapses  and  revivals  in  the 
echo  signal  owe  to  the  natural  abundance  (1.1%)  of  13 C  nuclei  in  the  lattice  environment. 
Free  to  rotate  in  its  lattice  site,  13C  experiences  Larmor  precession  due  to  the  static  magnetic 
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field  B0.  This  causes  a  periodic  fluctuation  in  the  interaction  between  the  nuclei  and  the 
NV  spin,  thus  leading  to  an  oscillating  phase  accumulation  during  the  echo  cycle  HI  71  j. 
If  the  free  precession  time  r/2  coincides  with  an  integer  multiple  of  the  Larmor  period  2%, 
then  the  net  phase  accumulation  from  the  13C  interaction  is  cancelled  out,  leading  to  a 
revival  in  the  fluorescence  signal.  Other  precession  time  intervals  lead  to  modulation  of 
the  echo  signal,  alternating  between  collapses  and  revivals.  Eventually,  weak  dipole-dipole 
interactions  between  pairs  of  13C  nuclei  and  between  13C  nuclei  and  spins  of  paramagnetic 
impurities  induce  flip-flops  in  the  interaction  between  the  NV  spin  and  13C  nuclei,  leading 
to  decoherence  |72|.  To  ensure  that  all  13  C  nuclei  precess  at  the  same  frequency,  special 
care  is  taken  to  ensure  zero  misalignment  of  the  static  magnetic  field  with  the  NV  symmetry 
axis.  A  transverse  static  field  induces  anisotropic-hyperfine  coupling  between  the  NV  spin 
and  the  13C  nuclei  |[73j,  causing  a  position-dependent  precession  frequency  and  thus  the 
disappearance  of  the  collapses  and  revivals,  and  a  reduction  in  the  decoherence  time  T2. 

The  envelope  of  the  signal  is  modelled  as  a  stretched  exponential  e~^TlT^v  with  stretching 
parameter  p.  For  a  single  NV  p  =  31m  but  for  ensembles  p  usually  varies  between  0.5  and 
3  [73].  The  curve  fits  for  our  diamonds  yielded  0.5  <  p  <  2.  The  collapses  and  revivals  in 
the  fluorescence  signal  due  to  the  13C  nuclei  is  theoretically  a  fourth-power  sinusoid  [75J  for 
a  single  NV  spin,  but  we  used  a  linear  combination  of  Lorentzian  and  Gaussian  shapes  for 
each  revival  to  modulate  the  stretched  exponential.  From  this  fit  we  extract  T2: 

s(t)  =  4,  +  A  exp[(— t/T2)p]  £  ( B9„e-(T-”''|2/(2G5)  +  - B‘n  )  ,  (4.2) 

ti  V  (r-  ntTf  +  (f )  ) 

where  Ao,  Ai,  Bgn,  Bln,  Gg,  Gi  are  fit  parameters,  tr  is  the  revival  time  of  first  revival,  and  n 
ranges  from  1  to  the  number  of  revivals  N  in  the  signal. 


4.4.2  Results 


Figure  T4|  gives  the  measured  T2  values  for  all  the  diamonds.  The  error  bars  are  extracted 
from  the  68%  confidence  intervals  obtained  from  the  decoherence  curve  fitting.  For  irradia¬ 
tion  doses  up  to  1  x  1018  e“/cm2,  T2  ranges  from  over  380  -  460  /rs;  after  this  dose  the  values 
begin  to  decline.  It  is  interesting  to  note  that  this  dose  is  also  where  the  florescence  and 
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Figure  4-4:  Spin-decoherence  times  T2  for  all  diamonds  as  a  function  of  irradiation  dose. 


charge  state  fluorescence  ratio  also  begin  to  change  behavior  as  seen  in  chapter  |Chapter~3l 
past  this  dose,  the  total  NV  fluorescence  begins  to  decline  and  the  charge  state  fluorescence 
ratio  begins  to  increase  relative  to  lower  doses  at  high  laser  excitation.  The  nitrogen  content 
(~  0.15  ppm)  and  13C  content  of  the  diamonds  is  the  same,  which  is  expected  to  be  the  reason 
for  the  uniformity  of  T2  at  low  doses.  Past  1  x  1018  e“/cm2,  however,  T2  is  no  longer  limited 
by  the  substitutional  nitrogen  or  13C  concentration  in  the  diamond.  Further  investigation 
is  needed  to  elucidate  the  nature  of  this  new  limiting  mechanism.  One  possible  cause  could 
be  the  creation  of  paramagnetic  vacancy-related  defects  at  these  high  doses  mm,  which 
couple  to  the  NV  center  to  induce  decoherence  at  a  faster  rate  than  14N  and  13C.  Electron 
paramagnetic  resonance  (EPR)  spectroscopy  is  a  technique  to  identify  paramagnetic  defects; 
measuring  EPR  spectra  for  the  high  dose  diamonds  and  comparing  to  the  literature  would 
indicate  what  defects  are  formed. 
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Figure  4-5:  (top)  Pulse  sequence  used  to  measure  the  spin-relaxation  time  T).  (bottom) 
Example  signal  from  a  spin-relaxation  pulse  sequence,  where  the  calculated  spin-relaxation 
time  is  T)  =  5.4  ±  0.0751  ms. 


4.5  Spin-Lattice  Relaxation  Time  T\  Measurements 


4.5.1  Measurement  Procedure 


To  measure  T)  time  scales,  the  pulse  sequence  depicted  in  |Figure  47q|  is  employed.  After 
being  pumped  to  the  |0)  state,  a  resonant  MW  ir  pulse  flips  the  spin  to  the  |  —  1)  state  and  is 
allowed  to  relax  toward  the  thermal  equilibrium  mixed  state  over  a  variable  time  r,  yielding 
the  fluorescence  trace  mi.  The  NV  ensemble  is  then  pumped  back  to  |0)  and  allowed  to 
relax  again  to  thermal  equilibrium  with  no  MW  pulse,  giving  the  fluorescence  trace  mo. 
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Figure  4-6:  Spin-lattice  relaxation  times  T\  as  a  function  of  dose.  All  diamonds  that  had 
Ti  measurements  under  high  laser  excitation  and  experienced  partial  re-polarization  during 
relaxation  are  omitted. 


4.5.2  Results  and  Limitations 


To  form  the  signal  and  calculate  spin-lattice  relaxation,  we  form  the  population  difference 
m0  —  mi .  This  follows  a  simple  exponential  profile  with  decay  time  T\.  Similar  to  the  spin- 
echo  and  Ramsey  measurements,  the  subtraction  of  the  signals  eliminates  common-mode 
noise  from  laser  fluctuations  during  the  measurement  and  the  common  fluorescence  from  the 
NV  orientations  not  addressed  by  the  MW  ir  pulse.  |Figure  T6|  shows  the  values  of  T\  as  a 
function  of  irradiation  dose.  The  error  bars  are  extracted  from  the  68%  confidence  intervals 
obtained  from  the  spin-relaxation  curve  fitting.  The  spread  in  error  bars  are  the  result  of  a 
spread  of  differing  measurement  times;  smaller  error  bars  correspond  to  more  averages  of  a  Ti 
measurement  (the  error  bars  themselves  arise  from  the  confidence  interval  exponential  fitting 
parameter  b  in  the  signal  s  =  ae~br).  No  clear  trend  can  be  seen  in  these  measurements, 
which  suggest  that  Ti  is  not  strongly  affected  by  irradiation  and  annealing. 
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4.6  Conclusion 


In  this  chapter,  we  have  described  the  role  that  various  decoherence  time  scales  play  in  NV 
magnetometry.  For  characterizing  the  diamonds  for  AC  magnetometry,  the  spin-decoherence 
time  T2  was  measured  with  a  spin-echo  pulse  sequence.  For  diamonds  in  this  study,  irra¬ 
diation  doses  greater  than  1  x  1018  e_/cm2  limit  coherence  beyond  the  limit  imposed  by 
substitutional  nitrogen  and  13C  alone.  Similarly,  for  DC  magnetometry,  we  measured  the 
spin-dephasing  time  T2*  using  a  Ramsey  pulse  sequence.  The  measured  T2*  times  fall  below 
the  1  /is  limit  imposed  by  13C,  either  due  to  systematic  effects  in  the  experiment  and/or  to 
effects  imposed  by  irradiation  and  annealing.  Further  measurement  of  Tf  and  T2*  will  be  the 
focus  of  future  work. 
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Chapter  5 


Analysis  and  Outlook 

5.1  Overview 

Using  the  fluorescence  and  coherence  measurements  discussed  in  previous  chapters  of  this 
thesis,  we  now  determine  the  optimum  diamond  processing  regime  for  commercially  available 
CVD  diamonds  suited  for  NV  ensemble  magnetometery.  A  simple  model  is  introduced  to 
calculate  the  magnetic  sensitivity  figure  of  merit  of  each  diamond.  A  discussion  of  the  results 
is  followed  by  proposed  follow-on  experiments  to  achieve  a  deeper  understanding  of  the  data. 


5.2  Sensitivity  Figure  of  Merit 

5.2.1  Modelling  the  Sensitivity 

So  far  in  this  thesis  we  have  only  discussed  the  behavior  of  the  fluorescence  and  coherence 
time  separately  as  a  function  of  optical  excitation  power  and  irradiation  dose.  The  original 
question  of  this  study,  however,  deals  with  their  combined  effect  on  magnetic  sensitivity  of 
an  NV  ensemble  magnetometer.  Recall  that  the  magnetic  sensitivity  equation  rj  l/a,\fWT 
depends  on  the  fluorescence  contrast  a,  photons  collected  during  measurement  /?,  and  the 
measurement  time  T  (optimized  when  set  on  the  order  of  T2*  f°r  DC  and  T2  for  AC  magne- 
tometry).  For  our  commercially  available  bulk  CVD,  high  purity  (150  ppb  [N5]),  electron 
irradiated  diamonds,  we  now  analyze  our  collected  data  to  see  which  irradiation  and  anneal- 
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ing  regime  optimizes  the  diamond  for  ensemble  NV  magnetometry. 

We  will  focus  only  on  the  main  variables  of  the  sensitivity  equation  r\  ~  1  /ay/JJT  and 
ignore  fundamental  constants  and  factors  of  order  unity  that  arise  from  the  specific  magne¬ 
tometry  scheme  implemented.  Recall  that  a  given  NV  ensemble  emits  a  photoluminescence 
(PL)  spectrum  under  excitation  with  contributions  from  both  the  NV°  and  NV-  charge 
states.  In  |Chapter~3]  we  modelled  these  contributions  as  arising  from  a  linear  combination 
of  the  individual  charge  state  PL  spectra,  f^v 0  and  Jnv-  ,  where  /  is  defined  at  the  total 
integrated  intensity  of  the  individual  PL  spectra.  In  this  chapter  we  take  these  values  to 
be  fluorescence  rates  in  photons/second.  The  total  fluorescence  f^v  is  then  taken  to  be 
a  linear  combination  of  the  two  rates:  fwv  =  co_/W°  +  C-Inv-i  where  cq  and  c_  are  the 


fitting  coefficients  from  the  linear  combination.  In  Chapter  1  the  contrast  a  was  defined  in 
terms  of  two  measured  signals  in  a  standard  NV  measurement,  m0  and  m i.  The  normalized 
difference  between  these  two  signals  provides  a  measure  of  the  population  difference  of  the 
NV  center  being  addressed.  For  a  single  NV-  center,  when  the  ground  state  spin  is  in  |  ±  1), 
the  maximum  population  difference  is  30  %  [7].  Assuming  this  occurs  under  the  conditions 
where  signal  mi  is  taken,  we  write  the  fluorescence  as  m0  =  fNV-  ,  mi  =  (1  —  0.3 )fNV-  R 
In  an  ensemble  measurement,  however,  the  presence  of  NV°  serves  to  decrease  the  contrast 
by  adding  background  fluorescence.  Furthermore,  as  explained  in  |Appendix  C[  typically  in 
NV  ensemble  magnetometry  only  one  NV  orientation  in  the  crystal  lattice  is  addressed  by 
microwave  (MW)  control  pulses,  while  the  other  three  classes  provide  additional  background 
fluorescence.  For  laser  polarization  normally  incident  on  a  diamond  with  (100)  face  geometry, 
only  two  NV  classes  can  be  preferentially  excited  at  a  time,  while  the  other  two  are  minimally 
excited.  For  this  geometry  it  has  been  measured  that  the  intensity  ratio  of  maximum  to 
minimum  fluorescence  is  2:1  176],  Therefore,  out  of  the  total  fluorescence  of  all  four  NV 
classes,  2/3  arises  from  the  three  NV  classes  (one  preferentially  excited  and  two  minimally 
excited)  not  addressed  and  1/3  arises  from  the  NV  class  (preferentially  excited)  that  is  being 
addressed.  When  the  latter  NV-  class  is  in  the  |  ±  1)  state,  the  fluorescence  from  the  NV- 
ensemble  is  tti\  =  Co/jvy°  T  c_  (|/atv-  V  (1  —  0.3)^/j\ry-)  =  Co/jvy°  T  0.9 c_/tvv-,  where  we 


1Note  that  30%  is  the  typical  maximum  contrast  that  can  be  achieved;  if  laser  broadening  exists  the 
contrast  will  decrease.  In  a  pulsed  experiment,  however,  this  is  not  an  issue  as  long  as  the  pulses  are  well 
calibrated  and  no  contrast  is  lost  to  decoherence. 


82 


have  included  the  NV°  fluorescence  that  is  the  same  for  m i  and  m0.  Putting  this  all  together, 
the  measurement  contrast  a  =  (mo  —  mi) /(mo  +  mi)  is 


_  O-lc-iW- _  /g 

2co/atv°  +  1-9  C-fNv- 

The  photons  collected  /3  is  modeled  as  m°2  ~tgate,  where  tgate  is  the  fluorescence  collection 
time.  In  terms  of  the  fluorescence  rates  calculate  above,  this  yields 


P  —  (cofNV0  +  0.95/W-)  tgate ■  (5-2) 

The  measurement  time  T  is  the  length  of  the  free-precession  time  in  the  measurement, 
which  we  take  here  to  be  T2,  and  the  dead-time  of  a  measurement  used  to  re-pump  the  NV 
ensemble  to  the  |0)  state.  The  sensitivity  formula  rj  ~  1  /a>/]3T  assumes  that  T2  tpump ;  if 
the  coherence  time  and  pump  time  are  comparable  then  the  sensitivity  needs  to  be  modified  to 
Tj  r^J  yr2  + 1„/  cry /73T2.  To  assist  in  comparing  relative  sensitivities  between  the  diamonds, 
we  use  the  inverse  sensitivity  as  our  figure  of  merit: 

FOM  =  a^T2  .  (5.3) 

\f  2  2  T  tpUrnp 

In  terms  of  the  figure  of  merit,  a  high  sensitivity  (low  value  of  rj)  corresponds  to  a  high 
FOM.  Thus,  we  now  seek  to  determine  which  processing  conditions  gives  the  highest  values 
of  FOM. 

To  calculate  FOM  for  each  diamond,  we  use  the  values  of  the  fitting  coefficients  c0,  C\ 
and  the  measured  spin-decoherence  time  T2  for  each  diamond  for  all  laser  powers  used.  One 
more  thing  is  needed,  however.  In  practice,  due  to  the  PL  spectra  overlap  of  NV-  with  NV°, 
a  long-pass  optical  filter  is  used  in  measuring  m0  and  mi  to  improve  the  contrast  a.  While 
this  does  prevent  some  NV-  fluorescence  from  being  measured,  the  much  larger  decrease 
in  NV°  fluorescence  compensates  for  this.  The  cut-off  wavelength  Ao  of  this  filter,  however, 
depends  on  the  relative  contributions  of  each  charge  state.  The  individual  charge  state 
total  fluorescence  f^v0  and  Jnv-  are  normalized  to  have  an  area  of  unity.  In  calculating 
the  FOM,  we  varied  the  cut-off  wavelength  Aq  for  all  wavelength  measured  (1340  with  our 
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Figure  5-1:  Sensitivity  figure  of  merit  FOM  as  a  function  of  long-pas  filter  wavelength  A0 
for  Rutherford  (1  x  1016  e“/cm2,  1250°C  anneal)  at  70  mW  of  excitation. 


CCD),  which  varied  the  value  of  the  total  fluorescence  f^v °-  from  0  to  1.  Figure  5-1  shows 
an  example  of  the  variation  of  FOM  as  a  function  of  the  cut-off  filter.  As  expected,  an 
optimum  value  is  reached  as  the  reduction  in  NV°  increases  the  contrast  a.  At  higher  cut-off 
wavelengths,  the  reduction  in  the  collected  photons  (3  from  both  charge  states  counteracts 
the  increase  in  contrast  and  reduces  the  FOM.  Out  of  the  1340  values  of  FOM  calculated, 
one  for  each  long-pass  filter,  we  report  only  the  maximum  values  in  the  rest  of  this  chapter. 


5.2.2  Sensitivity  vs.  Irradiation  Dose 


Figure  5^3|  plots  the  sensitivity  FOM  for  each  diamond  at  four  excitation  powers  for  both 


annealing  temperatures.  At  low  irradiation  doses  the  sensitivity  FOM  increases  until  about 
(~  1  x  1018  e"/cm2).  This  agrees  with  the  increasing  NV  fluorescence  seen  at  these  doses. 
The  drop  in  FOM  past  this  dose  agrees  with  the  drop  in  fluorescence  seen  in  |Figure  342]  and 
the  drop  in  T2  seen  in|Figure  4-4|  With  regards  to  the  optical  power,  the  FOM  monotonically 
increases  (ignoring  the  systematic  error  for  Bascom  at  (5  x  1017  e"/cm2,  850°C  anneal))  due 
to  the  increased  fluorescence  higher  excitation  brings.  The  relative  change  in  sensitivity, 
however,  decreases  as  the  power  increases,  which  agrees  with  the  start  of  saturation  seen  in 
Figure  3^6|  and  |Figure  3W[  as  well  with  the  decrease  in  contrast  due  to  a  lower  charge  state 
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Figure  5-2:  Optimum  long-pass  filter  wavelength  A0  as  a  function  of  laser  power  and  irradi¬ 
ation  dose  for  diamonds  annealed  at  (top)  850°C  and  (bottom)  1250°C. 
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Figure  5-3:  Relative  sensitivity  vs.  irradiation  dose  over  four  orders  of  magnitude  of  excita¬ 
tion  power  for  diamonds  that  underwent  (a)  850°C  annealing  treatment  (b)  1250°C  annealing 
treatment. 


fluorescence  ratio,  as  was  seen  in  Figure  3-9  (ignoring  the  highest  two  doses). 


With  regards  to  annealing  temperature,  jFigure  5-4|  shows  the  sensitivity  FOM  for  both 
anneal  temperatures  at  a  low  (100  /rW)  and  high  (150  mW)  laser  excitation.  At  low  laser 
power,  the  1250°  annealed  diamonds  achieve  the  optimum  FOM  in  the  dose  range  (1  x  1016 
e_/cm2  <  D  <  1  x  1018  e_/cm2),  which  is  consistent  with  the  total  fluorescence  plotted 
in  |Figure  3^3]  where  in  this  dose  range  1250°C  diamonds  have  a  higher  total  fluorescence 
than  850°C  annealed  diamonds.  At  higher  laser  powers,  diamonds  annealed  at  850°C  show 
higher  FOM,  in  agreement  with  the  lower  charge  state  fluorescence  ratio  for  1250°C  annealed 
diamonds  seen  in|Figure  3-10, 


5.3  Outlook 

The  data  presented  in  the  previous  section  reveal  that  for  the  majority  of  laser  powers  used, 
diamonds  annealed  at  1250°C  and  irradiated  with  a  dose  around  5  x  1017  e“/cm2  would 
give  an  NV  ensemble  magnetometer  the  highest  sensitivity.  At  all  laser  powers,  the  dose 
range  (5  x  1016  e~/cm2  <  D  <  5  x  1018  e_/cm2)  yields  the  optimum  sensitivity  compared 
to  other  doses  at  1250°C.  At  higher  laser  powers,  diamonds  annealed  at  850°C  with  a  dose 
around  1  x  1016  e“/cm2  optimize  the  sensitivity.  The  highest  sensitivity  FOM  obtained 
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Figure  5-4:  Relative  sensitivity  vs.  irradiation  dose  at  all  anneal  temperatures  for  (left)  100 
[i W  and  (right)  150  mW  of  laser  excitation. 


is  at  200  mW  for  850°C  annealing  temperature  and  a  irradiation  dose  of  1  x  1016  e_/cm2. 
Although  5  x  1018  e“/cm2  gives  a  slightly  higher  sensitivity,  the  increase  in  the  magnitude 
of  the  charge  state  ratio  seen  in  |Figure  3-1 1|  as  well  as  the  lower  value  of  T2  relative  to  the 
low  dose  diamonds  seen  in  |Figure  4A|  make  it  suspect  to  the  unknown  effect  that  the  highly 
irradiated  diamonds  have.  At  lower  powers  this  dose  not  yield  highest  FOM.  The  dose 
range  (1  x  1016  e“/cm2  <  D  <  1  x  1017  e“/cm2),  however,  consistently  gives  the  optimum 
sensitivity  FOM  regardless  of  laser  power. 

This  answers  the  question  this  thesis  originally  set  out  to  answer.  However,  as  seen  in 
the  data,  at  high  doses  (>  1  x  1019  e"/cm2)  unobserved  deviations  are  seen  from  expected 
behavior  in  terms  of  a  simple  model  of  monovacancies  and  substitutional  nitrogen  atoms 
in  the  diamond  lattice.  At  these  higher  doses  the  radiation  damage  to  the  lattice  creates 
additional  defects  that  are  speculated  to  degrade  the  spin  decoherence  time  T2  of  the  diamond 
and  significantly  affect  the  charge-state  dynamics  of  the  NV  center,  shown  by  the  increase 
in  charge  state  fluorescence  ratio  under  increased  excitation  power.  If  the  results  of  this 
study  are  to  be  generalized  to  a  wider  range  of  diamonds,  to  include  differing  nitrogen 
concentrations,  it  is  important  that  the  behavior  seen  at  high  dose  is  fully  understood. 

The  behavior  of  the  charge  state  properties  and  coherence  times  under  high  dosages  is 
most  likely  due  to  the  presence  of  increased  vacancy- related  defects  [38,  55J,  such  as  vacancy 
chains.  Since  a  majority  of  these  defects  are  paramagnetic,  they  can  be  identified  by  electron 
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Figure  5-5:  Example  EPR  spectrum  for  Ladd-Franklin  (5  x  1018  e“/cm2,  850°C  anneal). 
While  the  signal  ideally  should  have  a  flat  baseline,  the  MW  power  causes  a  response  from 
holding  tube  that  needs  to  be  subtracted. 


paramagnetic  resonance  spectroscopy.  Briefly,  under  a  fixed  microwave  frequency  (around 
10  GHz  for  typical  X-band  EPR)  an  applied  static  magnetic  field  is  swept  through  a  range  of 
amplitudes,  altering  the  Zeeman  energy  splittings  of  unpaired  spins  in  paramagnetic  defects. 
When  the  magnitude  of  the  field  causes  the  splitting  to  match  the  applied  MW  frequency,  a 
transition  between  spin  states  occurs  and  the  absorption  of  the  MW  is  recorded.  To  enhance 
the  signal-to-noise,  phase  sensitive  detection  is  used  where  the  magnetic  field  is  modulated 
as  it  is  swept.  This  results  in  the  measurement  of  the  first  derivative  of  the  actual  MW 
absorption  signal.  |Figure  5^5|  shows  an  example  raw  EPR  spectra  for  a  diamond  irradiated 
with  a  dose  of  5  x  1018  e~/cm2.  The  uneven  baseline  is  the  EPR  signal  from  the  Roxolite 
tube  used  to  hold  the  diamond  in  the  MW  cavity.  Comparison  with  Figure  2  in  [38j  indicates 
that  some  of  the  defects  created  at  this  high  dose  are  vacancy  chains.  A  systematic  study 
of  EPR  spectra  for  the  diamonds  in  this  thesis  could  identify  paramagnetic  defects  created 


at  high  electron  irradiation  doses,  which  could  assist  in  understanding  the  behavior  of  the 
charge  state  fluorescence  and  coherence  times  at  these  high  doses. 
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Appendix  A 


Experimental  Setup  and  Calculations 


This  appendix  describes  the  experimental  setup  used  for  this  thesis  and  provides  technical 
details  relevant  to  the  acquisition/analysis  procedures  discussed  in  earlier  chapters.  In  ad¬ 
dition,  the  laser  intensity  for  a  given  measured  laser  power,  the  approximate  region  of  the 
confocal  spot  that  is  imaged  by  the  collection  fiber,  and  the  approximate  laser  power  needed 
to  reach  NV-  saturation  is  calculated. 


A.l  Experimental  Setup 

A.  1.1  Overview  of  the  Setup 

A  home-built  fluorescence  microscope  served  as  the  basis  for  both  photoluminescence  (PL) 
spectroscopy  and  coherence  time  measurements  of  the  diamond  samples.  The  apparatus  is 
shown  in |Figure  AT. 
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Figure  A-l:  Diamond  characterization  experimental  setup  (top)  photograph  and  (bottom) 
schematic. 
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The  two  main  sections  of  this  setup  are  the  optical  excitation  path  and  the  NV  fluores¬ 
cence  collection  path.  The  optical  excitation  path  begins  with  a  5  W  532  nm  laser  (Coherent 
Verdi  G5).  To  maximize  the  proportion  of  NVs  in  the  negative  charge  state  the  optimal  ex¬ 
citation  wavelength  is  in  the  range  510  -  540  nm  [19J,  even  though  the  NV-  zero-phonon 
line  resonance  is  at  637  nm.  To  control  the  intensity  of  light  from  the  laser  box  into  the 
setup  the  green  laser  light  first  traverses  a  half-wave  plate  (HWP)  and  a  polarizing  beam 
splitter  (PBS).  Amplitude-controlled  light  is  deflected  by  an  acousto-optic  modulator  (AOM; 
Crystal  Tech/Gooch  &  Housego  3200-147,  driver  model  1200-200).  The  first-order  diffracted 
beam  passes  through  a  second  HWP  and  PBS,  which  is  used  to  control  the  power  into  a 
microscope  objective.  The  beam  then  passes  through  a  third  HWP  that  establishes  reference 
polarization  axes  of  the  light  into  the  objective.  The  beam  reflects  off  a  dichroic  mirror  (Sem- 
rock  LPD02-532RU,  edge  wavelength  537.2  nm)  to  a  microscope  objective  (Mitutoyo,  20x, 
NA  =  0.42,  infinity-corrected)  and  is  focused  to  a  spot  size  3.01  /im  in  diameter  (calculated 
below).  The  diamond  is  located  on  a  ^-translation  stage  directly  over  the  objective  shown  in 


Figure  A-2  (left).  A  loop  antenna  connected  to  an  amplifier  (Mini-Circuits  ZHL-16W-43+) 


and  vector  signal  generator  (Keysight  N5182B)  provides  microwaves  for  coherent  control  of 
the  NV  ensemble.  A  1"  pedestal  spacer  serves  as  the  diamond  mount;  glass  coverslip  slides 
glued  onto  the  spacer  hold  the  diamond  in  a  consistent  location  while  under  pressure  from 
the  loop  antenna  as  shown  in  Figure  A-2  (right)]  For  coherence  measurements,  the  diamond 
mount  is  aligned  such  that  the  confocal  spot  is  at  the  top  of  the  diamond  near  the  antenna 
to  provide  a  uniform  MW  field,  while  for  PL  measurements  the  mount  is  adjusted  such  that 
the  confocal  spot  is  at  the  bottom  of  the  diamond  closest  to  the  objective  to  avoid  laser  or 
NV  fluorescence  absorption.  To  provide  a  static  bias  magnetic  field,  a  permanent  magnet  is 
attached  to  a  translation  stage  behind  the  diamond  mount. 


Light  emitted  from  the  NV  ensemble  in  the  diamond  is  collected  by  the  objective  and 
is  transmitted  through  the  dichroic  mirror  to  a  0.10  fi m  core  multi-mode  fiber  (Thorlabs, 
M64L02).  This  coupling  element  is  attached  to  a  steel  translation  stage  and  serves  as  a 
pinhole  to  isolate  the  NV  fluorescence  emitted  from  the  confocal  volume  of  the  objective  in 
the  diamond.  This  fiber  is  connected  to  one  of  two  other  fibers,  a  spectrometer  for  photolu¬ 
minescence  measurements  or  a  single  photon  counter  for  coherence-time  measurements. 
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Figure  A-2:  Diamond  mount,  (left)  Vertical  view  of  the  objective-diamond  assembly.  Laser 
light  reflects  up  to  the  objective  in  a  cage  mount  where  it  focuses  onto  the  diamond  mount 
translation  stage,  which  is  adjustable  to  ensure  that  the  confocal  volume  is  at  the  top  of  the 
diamond  near  the  MW  loop  for  coherence  measurements  or  near  the  bottom  of  the  diamond 
for  PL  spectroscopy  measurements,  (right)  Left:  photograph  of  the  diamond  mount.  The 
loop  antenna  is  pressed  against  the  face  of  the  diamond  and  centered  about  the  confocal  spot 
to  maximize  MW  uniformity.  Right:  diamond  mount  schematic  with  top  and  side  views. 


A. 2  Fluorescence  Calculations 


In  this  section  we  calculate  the  relevant  experimental  parameters  for  the  setup  defined  above. 
We  will  model  the  intensity  of  the  green  laser  beam  focused  by  the  objective  as  a  Gaussian 
profile: 

2  P  f  2 r2  \ 

I(r,z)  = - rw^expl - —  ),  (A.l) 

where  P  is  the  laser  power,  w(z)  =  Woy/l  +  ( ~z]~zr )2  is  the  1/e2  beam  radius  (defined  as  the 
distance  where  the  intensity  drops  to  13.5%  of  its  maximum  value),  zr  =  irridwl/X  is  the 
Rayleigh  range  (adjusted  for  the  diamond  index  of  refraction  rid),  r  is  the  radial  distance  in 
the  plane  transverse  to  the  beam,  and  z  is  the  axial  distance  with  the  origin  defined  at  the 
center  of  the  confocal  volume  [77J.  The  beam  waist  diameter  at  the  origin  is  2 wq.  The  peak 
intensity  %  is  defined  to  be  Iq  = 
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Figure  A-3:  Profile  of  the  laser  beam  right  before  being  reflected  off  the  dichroic.  The  D4<rX 
and  D4ctY  give  measures  of  the  beam  widths  in  the  X  and  Y  directions. 


A. 2.1  Power  Density 


To  calculate  the  intensity  of  the  light  in  the  confocal  volume  of  the  diamond,  we  first  measure 
the  spot  size  of  the  green  laser  beam  before  the  dichroic  with  a  beam  profiler  (Newport  LBP2- 


HR-VIS).  Figure  A-3  shows  the  values  of  D4aX  and  D4<rY  of  the  beam,  which  for  a  Gaussian 
beam  give  the  width  in  terms  of  the  1/e2  distanc^j]  The  Verdi  G5  laser  manual  specifies  a 
Gaussian  beam  diameter  of  2.25  mm  ±  10%,  which  is  in  agreement  with  the  X-width  of  2.33 
mm  and  Y  width  of  2.24  mm.  Possible  reasons  for  the  different  X  and  Y  lengths  are  the 
diffraction  of  the  beam  from  the  AOM  and  imperfect  optics. 

From  the  Mitutoyo  objective  manual,  for  an  incident  Gaussian  beam  of  (1/e2)  width 
Diaser,  the  focused  beam  waist  is  Dsvot  =  For  an  input  radius  of  2.25  mm,  laser 

^  'R-LJla - 


wavelength  A  =  532  nm,  and  objective  focal  length  /0y  =  10  mm,  the  beam  waist  is  3.01  /im. 
This  yields  a  Rayleigh  range  zr  =  32  ym  for  a  confocal  volume  of  about  450  /im3.  Incident 
laser  power  was  measured  prior  to  the  dichroic  mirror.  Accounting  for  83%  transmission 


1D4(jX  is  defined  as  four  times  the  second  x-moment  Mx  of  the  transverse  intensity  distribution  of  the 
light  beam.  The  second  moment  in  the  x  direction  is  Mx 


^  f  f  I(x,y,z)dxdy 


T7  r2  cos2  0I(r,z)rdrd0  _ 

f  f  I (r,z)rdrd6 

w(z)/2.  Since  w  is  the  1/e 2  radius  and  D4crX=  AMX  =  2 re,  for  a  Gaussian  beam  both  D4crX  and  the  1/e2 
diameter  are  the  same. 
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of  the  dichroic  and  17%  front  surface  reflection  at  normal  incidence  to  the  diamoncQ  the 
intensity  of  the  laser  at  the  middle  of  the  confocal  spot  in  the  diamond  is 


(A.2) 


where  the  measured  laser  power  prior  to  the  dichroic  mirror  P  is  in  mW. 

A. 2. 2  Saturation  Power 

The  saturation  intensity  of  the  NV-  center  is  approximately  100  kW/cm2  |78|.  Thus,  for  a 
confocal  spot  size  diameter  of  3.01  pm  and  69%  power  loss,  this  corresponds  to  a  measured 
laser  power  of  5  mW  before  the  dichroic  in  this  experiment. 

A. 2. 3  Imaging  Region  Calculations 

In  |Chapter~3|  saturation  of  the  NV  fluorescence  is  considered.  To  see  saturation  effects  in 
confocal  microscopy,  it  is  important  that  the  region  imaged  by  the  pinhole  is  restricted  to  the 
central  region  of  the  confocal  spot,  where  saturation  is  noticeable  (i.e.  to  avoid  measuring 
the  increasing  fluorescence  in  the  wings  of  the  NV  fluorescence  profile  which  do  not  saturate 
as  fast).  The  collected  fluorescence  from  the  objective  used  the  full  numerical  aperture. 
With  a  back  aperture  of  8.4  mm,  the  collection  region  of  fluorescence  by  the  objective  is 
D collection  =  n  4A^o6j - ,  which  for  a  wavelength  range  of  560  -  818  nm  is  0.85  to  1.24  pm. 

Tf  Ubackaperture 


2 For  light  at  normal  incidence  to  a  diamond-air  interface,  the  reflection  is 


nd-nair 


Fresnel  equations,  where  nair  =  1.000  is  the  index  of  refraction  for  air  and  rid  =  2.419  is  the  index  of 
refraction  for  diamond. 
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Appendix  B 


NV  Fluorescence  Methods  and 
Calculations 


Chapter  ~3|  showed  the  results  of  fluorescence  measurements  taken  on  the  diamonds  as  the 
irradiation  dose  and  optical  excitation  power  were  varied.  To  extract  the  charge  state  flu¬ 
orescence  ratio  and  the  total  fluorescence  from  each  charge  state,  the  fluorescence  fraction 
of  a  photoluminescence  (PL)  spectrum  from  each  charge  state  must  be  calculated.  This  ap¬ 
pendix  details  the  charge-state  fluorescence  data  acquisition  process  and  the  accompanying 
analysis  to  calculate  the  charge  state  fluorescence  ratio  and  total  fluorescence  from  a  linear 
combination  of  NY0  and  NV-  PL  spectra. 


B.l  Charge  State  Data  Acquisition 

B.l.l  Measurement  Procedure 


All  PL  spectroscopy  measurements  are  performed  with  a  spectrometer  (Princeton  Instru¬ 
ments,  IsoPlane  SCT  320)  and  an  attached  CCD  camera  (Princeton  Instruments,  PIXIS 
100),  shown  in  Figure  B-la|  As  explained  in  Appendix  C[  the  half-wave  plate  before  the 
dichroic  in  the  experimental  setup  is  adjusted  so  that  the  532  nm  excitation  equally  ex¬ 
cites  all  four  NV  orientation  classes.  Princeton  Instruments  Intelical  USB-powered  sources 
provide  wavelength  and  intensity  calibrations  to  the  recorded  spectra.  The  intensity  calibra- 
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(a)  (b) 


Figure  B-l:  Experimental  apparatus  for  PL  spectroscopy  (a)  Spectrometer.  A  fiber  is 
mounted  to  the  entrance  slit  where  the  acquired  NV  fluorescence  is  diffracted  and  pro¬ 
cessed  by  the  CCD  camera.  The  position  of  the  fiber  can  be  adjusted  in  two  dimensions. 
A  micrometer  allows  fine  adjustment  of  the  spectrometer  input  slit,  (b)  Fiber  shaker  to 
eliminate  modal  noise  in  the  acquired  PL  spectra. 


tion  removes  instrumental  artifacts  in  the  PL  spectra  that  arise  from  the  differing  spectral 
responses  from  various  optical  elements  in  the  spectrometer  (mirrors,  diffraction  grating, 


CCD  camera,  etc.).  Figure  B-2  (top)  shows  a  sample  PL  spectrum  with  and  without  this 
calibration.  Only  the  rows  of  the  CCD  illuminated  by  the  fiber  are  used  in  calibrating  the 
spectra.  To  obtain  maximal  signal-to-noise,  the  exposure  time  of  each  frame  is  set  so  that 
the  peak  intensity  on  the  CCD  is  about  55,000  counts;  this  ensures  that  the  CCD  operates 
in  the  linear  regime  but  that  the  pixel  wells  are  still  nearly  full.  The  diffraction  grating  is 
set  to  300  gratings/mm  with  a  center  wavelength  of  690  nm,  which  records  spectra  in  the 
range  of  560  -  818  nm.  Finally,  to  avoid  signal  to  noise  degradation  by  modal  noise  in  the 
multi-mode  fiber  pi,  the  fiber  is  taped  to  a  speaker-head  driven  by  a  200  Hz  sinusoidal 


voltage  (15  Volts  peak-to-peak).  Figure  B-2  (bottom)  shows  a  sample  PL  spectrum  with 
and  without  the  fiber  shaker,  revealing  the  enhancement  in  spectral  resolution  of  the  NV- 
zero-phonon  line  amidst  the  NV°  phonon  sideband. 


After  data  acquisition,  the  resulting  files  are  exported  to  .csv  files  to  be  read  by  MATLAB 
code.  All  frames  taken  for  a  single  measurement  are  averaged  together,  followed  by  normal¬ 
ization  by  the  exposure  time  of  the  measurement.  To  analyze  the  charge  state  properties  of 
the  PL  spectra  for  each  diamond,  we  employ  a  linear  combination  of  NV°  and  NV-  spectra 
|65j  as  described  in  the  next  section. 


98 


Figure  B-2:  Spectrum  modification  to  correct  for  experimental  artifacts,  (top)  Normalized 
photoluminescence  spectra  of  the  diamond  Hilbert  (5  x  1015  e_/cm2,  1250°C  anneal)  at 
1  mW  optical  excitation  with  and  without  the  intensity  calibration  turned  on.  (bottom) 
Photoluminescence  spectra  of  the  diamond  Anderson  (1.6  x  1019  e“/cm2,  1250°C  anneal) 
at  300  //W  optical  excitation  with  the  fiber  shaker  turned  off  and  on.  Notice  how  the  NV- 
zero-phonon  line  is  resolvable  with  the  fiber  shaker  on. 
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B.1.2  Linear  Combination  of  NV  Spectra 

The  other  method  employed  to  analyze  the  charge  state  ratio  was  fitting  each  acquired  PL 
spectra  to  a  linear  combination  of  NV-  and  NV°  ‘basis’  spectrum.  That  is,  for  a  normal¬ 
ized  NV-  basis  function  I_( A)  and  normalized  NV°  basis  function  /o(A),  each  acquired  PL 
spectrum  I  is  fit  to  a  linear  combination  of  the  basis  functions 

/(A)  =  c_/_(A)  +  c0/0(A),  (B.l) 

where  c_  and  c0  are  the  weights  of  each  basis  function  whose  ratio  yields  the  desired  charge 
state  ratio. 


Constructing  the  Basis  Functions 

To  acquire  the  NV°  basis  function,  we  first  measured  the  PL  spectrum  of  Bascom  (5  x  fO17 
e“/cm2,  850°C  anneal)  at  200  mW  excitation  power  (the  highest  power  we  can  obtain  before 
exceeding  the  damage  threshold  of  the  AOM  in  our  setup).  This  PL  spectrum,  while  not 
completely  devoid  of  a  NV-  PL  fingerprint,  is  the  highest  fluorescing  diamond  in  our  diamond 
set  and  has  little  if  any  trace  of  a  NV-  ZPL,  so  to  a  good  approximation  can  serve  as  our 
NV°  basis  function.  To  obtain  the  NV-  basis  function,  we  measured  the  PL  spectrum  of 
Heaviside  (unprocessed)  at  fOO  //W.  The  low  power  ensures  that  this  specific  PL  spectrum 
is  predominantly  composed  of  NV-  fluorescence.  However,  since  this  spectrum  does  have  a 
NV°  PL  signature,  we  scale  the  NV°  basis  function  and  subtract  it  off  Heaviside’s  spectrum 


to  give  the  NV-  basis  function,  shown  in|Figure  B-3[  Both  basis  functions  are  then  divided 
by  their  integrated  area. 

To  check  the  validity  of  using  the  Bascom  PL  spectrum  as  the  NV°  basis  function,  we 
temporarily  mounted  a  405  nm  laser  pointer  and  aligned  it  to  excite  the  Bascom  diamond 
through  the  objective;  the  measured  spectrum  is  plotted  with  the  Bascom  spectrums  used 


for  the  NV°  basis  function  (excited  at  the  usual  532  nm)  in  Figure  B-4,  At  wavelengths 
below  490  nm  the  NV-  center  can  be  directly  ionized  to  NV°  with  a  single  photon  (as 
opposed  to  the  two-photon  process  outlined  in  |Chapter~H  which  results  in  a  lower  steady 
state  probability  to  be  in  the  negative  NV  charge  state  [19,  i2T].  To  within  a  small  error  the 
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Figure  B-3:  (Left)  PL  spectra  for  Bascom  (5  x  1017  e“/cm2,  850°C  anneal)  at  200  mW 
(top)  and  Heaviside  (unprocessed)  at  100  /iW  (bottom).  Bascom  serves  as  the  NV°  basis 
function  due  to  its  high  fluorescence  and  negligible  NV-  content.  (Right)  Construction  of 
the  NV-  basis  function.  The  NV°  basis  function  is  scaled  to  fit  the  NV°  ZPL  of  Heaviside’s 
PL  spectrum,  which  is  then  subtracted  to  yield  the  NV-  basis  function.  The  original  and 
basis-fitted  Heaviside  spectra  are  also  shown. 


two  functions  have  a  similar  spectral  profile. 


Calculating  the  Coefficients 

To  obtain  the  coefficients  c_  and  Co,  we  use  the  method  of  least-squares  regression  [80j. 
To  fit  the  measured  PL  spectrum  /(A)  to  the  form  /(A)  =  c_J_( A)  +  coio(A)  we  minimize 
the  sum  of  the  squares  of  the  residuals  JT(/(Aj)  —  c_/_(A*)  —  c0/o(Aj))2,  where  A j  are  the 
discrete  wavelengths  the  spectrum  is  measured  at.  This  yields  the  following  relations  for  the 
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Wavelength  (nm) 


Figure  B-4:  PL  spectra  for  Bascom  (5  x  1017  e  /cm2,  850°C  anneal)  at  200  mW  for  405  nm 
and  532  nm  excitation. 


102 


coefficients: 


YJi-YdhY-YJ-hYJh 
Y ,(i-YY,(hY-(Y,i-hY  ' 

YJ-YY AhY  -  (YJ-hY 


where  for  clarity  the  functional  dependence  on  \  is  suppressed.  The  error  of  this  fit,  which 
is  used  to  determine  the  error  of  the  charge  state  fluorescence  ratio  and  total  fluorescence, 
is  obtained  from  the  following  equations: 


<?i 


<t_  —  an 


(Jo  —  <*i\ 


ArraE^-c-7--^)2’ 

i 

E^o)2 

E,(/-)2E,(/o)2-(E^-/o)2’ 

E,(/-)2 

E?(/-)2E,(/o)2-(E?/-/o  )2’ 


where  N  is  number  of  wavelengths  sampled  (N  =  1340  for  this  CCD). 


(B.4) 

(B.5) 

(B.6) 


Ratio  and  Fluorescence  Calculations 


To  calculate  the  charge  state  fluorescence  ratio  R,  we  take  the  ratio  of  the  calculated  coeffi¬ 
cients  from  the  previous  section: 

R=—>  (B.7) 

c0 

where  the  error  is  given  by  standard  error  analysis: 


Sr  =  R 


/ 


2 

+ 


(B.8) 


The  fluorescence  of  each  charge  state  is  simply  the  coefficient  for  each  basis  function 
because  the  basis  functions  have  an  area  of  unity:  /_  =  c_,/0  =  c0,fl_  =  a_,fl0  =  cr0.  The 
total  fluorescence  is  then  given  by  the  sum  of  the  coefficients: 


iW  =  /o  +  /- , 


(B.9) 
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where  the  error  is  the  sum  of  the  coefficient  errors 


Sf  — 


(BIO) 
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Appendix  C 


Coherence  Time  Measurement 
Procedure 


Chapter  ~4|  presented  the  results  of  coherence  time  measurements  for  T|,T2,  and  Tf.  Before 


each  of  these  measurements  are  taken,  however,  certain  experimental  procedures  are  exe¬ 
cuted  to  ensure  the  measured  coherence  times  are  limited  by  the  concentration  of  nitrogen 
and  other  impurities  in  the  diamond  lattice  and  not  by  experimentally  controllable  sources 
of  decoherence.  These  include  anisotropic-hyperhne  coupling  to  the  13C  nucleus  by  the  NV 
center  due  to  a  misaligned  static  field  and  poorly  calibrated  proper  pulse  lengths  and  fre¬ 
quencies.  In  this  appendix  we  explain  the  experimental  methods  used  before  coherence  time 
characterization.  Specifically,  the  methodology  and  experimental  acquisition  of  continuous- 
wave  magnetically  detected  optical  resonance  (CAW-ODMR),  Rabi  nutation,  and  re-pump 
optimization  measurements  will  be  explained. 


C.l  Electronics 

For  coherence  time  measurements,  the  NV  fluorescence  collected  by  a  10-/rm  multimode  fiber 
is  sent  to  an  avalanche  photodiode  single-photon  counter  (Excelitas  SPCM-AQRH-FC).  The 
registered  counts  are  recorded  by  a  DAQ  (National  Instruments,  BNC-2120),  which  is  in 
turn  gated  by  a  pulse  generator  (SpinCore  Pulse  Blaster  ESR-PRO-500).  The  Pulse  Blaster 
controls  the  timing  of  the  optical  excitation  (by  means  of  switching  the  AOM  on/off),  the 
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Figure  C-l:  Wiring  diagram  for  coherent  control  of  the  NV  center.  PBO  controls  collection 
of  photons;  PB2  controls  the  MW  field  delivered  to  the  diamond;  and  PB3  switches  the  laser 
to  excite  and  polarize  the  NV  center. 


MW  excitation  through  the  loop  antenna,  and  the  measurement  of  the  photon  counts  via 
the  DAQ.  [Figure  C-l|  shows  the  wiring  of  the  SPCM,  AOM,  and  MW  drive;  the  “PB”  labels 
correspond  to  different  channels  of  the  Pulse  Blaster  card  controlling  each  instrument. 


C.2  Static  Magnetic  Field  Alignment 

C.2.1  Initial  Considerations 

The  crystal  structure  of  diamond  allows  the  NV  center  to  be  aligned  along  one  of  four  crystal¬ 
lographic  axes,  giving  rise  to  four  ‘classes’,  or  orientations,  of  NV  centers.  In  magnetometry 
schemes,  the  measurement  sequences  employed  usually  focus  on  only  one  class  of  NV  cen¬ 
ters.  To  isolate  control  of  this  NV  class  from  the  other  three  orientations,  a  static  magnetic 
field  is  applied  to  break  the  degeneracy  of  the  |  ±  1}  spin  manifold,  where  the  splitting  is 
proportional  to  the  projection  of  the  field  along  a  given  NV’s  axis  to  zeroth  order.  Due  to  the 
natural  abundance  of  13C  in  our  samples  (1.1%),  the  magnetic  field  must  be  aligned  to  the 
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desired  NV  class.  If  a  component  of  the  static  field  is  transverse  to  this  NV  axis,  anisotropic- 
hyperfine  coupling  between  the  NV  electronic  spin  and  proximal  13C  nuclei  degrades  the 
ensemble  coherence  time  m  The  degree  of  this  ensemble  coherence  time  degradation  is 
dependent  on  the  strength  of  the  transverse  field.  Thus,  stronger  applied  bias  fields  need  to 
be  better  aligned  to  the  NV  axis  compared  to  lower  applied  bias  fields. 

The  field  magnitude  along  the  NV  axis  is  another  important  consideration.  It  must  not 
be  too  low:  Ref.  l»l]  illustrates  that  below  10  G,  the  contribution  to  NV  decoherence  from 
13C  nuclei  increases,  thus  degrading  the  value  of  T2  for  the  diamond.  Most  importantly, 
though,  is  the  effect  it  has  on  the  signal  we  expect  from  a  spin-echo  measurement.  The 
Larmor  precession  of  the  13C  nuclei  cause  collapses  and  revivals  in  the  signal  from  a  spin- 
echo  measurement  as  discussed  in  Chapter  ~4|  The  period  of  the  Larmor  precession  (and 
thus  the  period  of  the  observed  collapses  and  revivals)  is  Tl,i3C  =  l/(7i3C-Bo),  where  B0  is 
the  magnitude  of  the  field  experienced  by  the  NV  ensemble  and  713c  is  the  gyromagnetic 
ratio  of  the  13C  nucleus  (1.071  kHz/G).  In  a  spin-echo  measurement,  the  fluorescent  signal 
experiences  a  revival  when  the  free  precession  time  r/2  between  the  (MW)  ir/2  and  n  pulses 
is  equal  to  an  integer  number  of  the  13C  Larmor  period.  To  have  enough  points  for  an 
accurate  fit  of  the  signal  we  chose  the  field  to  induce  at  least  10  revivals  in  the  signal,  which 
means  the  total  pulse  sequence  measurement  time  r  must  be  equal  to  T2/10.  Given  typical 
T2  for  diamonds  with  [ N ]  ~  0.15  ppm  and  13C  ~  1.1%  are  around  400  (as,  this  corresponds  to 
about  50  G.  We  have  assumed  here  that  the  field  doesn’t  vary  considerably  over  the  confocal 
volume  so  that  all  13C  nuclei  (barring  any  anisotropic  hyperfine  coupling  due  to  transverse 
fields)  precess  at  the  same  rate. 


C.2.2  (100)  Diamond  Geometry  and  Initial  Field  Alignment 

To  begin  field  alignment,  we  place  a  diamond  on  the  mount  such  that  the  magnet  is  roughly 
aligned  along  an  NV  axis,  as  shown  in  |Figure  Cb2|  This  arrangement  is  preferred  due  to 
the  (100)  face  and  sides  geometry  of  the  diamond.  The  MW  loop  is  positioned  on  top 
of  the  diamond  with  the  focused  laser  spot  in  the  center  of  the  antenna  loop.  The  loop  is 
pressed  against  the  diamond  to  maximize  the  MW  power  delivered  to  the  NV  ensemble.  The 
diamond  mount  is  adjusted  to  focus  the  confocal  spot  near  the  top  of  the  diamond  where 
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Green  =  vacancy  pair  1 
Blue  =  vacancy  pair  2 


Figure  C-2:  Diagram  of  the  dipole  magnet  and  (100)  diamond  geometry.  (Left)  Top  view 
of  the  magnet  and  diamond.  The  diamond  is  mounted  in  the  setup  such  that  the  corner 
is  bisected  by  the  dipole  axis.  (Right)  Side  view  of  the  diamond  and  magnet.  The  ideal 
alignment  for  this  setup  and  diamond  orientation  is  to  have  the  held  axis  aligned  along  the 
NY  class  shown. 


the  MW  loops  is  located  (see  Figure  A-2  (left)). 


This  arrangement  is  to  optically  address 


and  detect  only  the  NVs  that  are  near  the  antenna  loop  to  have  strong  MW  drive  held. 


C.2.3  CAW-ODMR  Measurements 


Continuous-wave  optically  detected  magnetic  resonance  (CAW-ODMR)  measurements  are 
performed  to  extract  the  resonances  of  the  NV  ensemble.  As  explained  in  |Chapter~Tj  in 
CAW-ODMR,  while  under  continuous  laser  excitation,  the  frequency  of  an  applied  MW  held 
is  swept  through  a  range  containing  the  NV  ensemble  resonances.  When  the  MW  held  is  not 
on  an  NV  resonance,  the  huorescence  intensity  is  a  maximum  as  the  NV  is  polarized  in  the  |0) 
state.  If  the  MW  is  on  an  NV  resonance,  it  drives  the  NV  population  which  cycles  between  the 
|0)  and  |±1)  states.  This  causes  a  reduction  in  the  huorescence  intensity  due  to  the  increased 
probability  of  these  states  to  decay  using  the  non-radiative  intersystem  crossing.  Thus,  the 


recorded  huorescence  signal  exhibits  dips  at  the  NV  resonance  frequencies.  |Figure  CV3 
(top) |  shows  the  pulse  sequence  used  to  obtain  CAW-ODMR  data.  For  the  hrst  huorescence 
signal  mi,  the  laser  and  MW  are  applied  continuously  while  the  “gate”  pulse  from  the  SPCM 
acquires  a  huorescence  signal.  In  the  second  half  of  the  sequence,  the  MW  is  switched  off  and 
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the  fluorescence  m0  is  recorded.  This  serves  as  the  baseline  which  the  signal  mi  is  normalized 
against,  forming  the  CAW-ODMR  contrast  mi /mo.  In  practice,  we  randomize  the  order  of 
the  frequencies  applied  to  eliminate  systematic  noise  (laser  intensity  drift,  thermal  drift, 
etc.). 

If  no  static  magnetic  field  is  applied,  the  only  transition  detected  is  the  degenerate  |0)  B 
|  ±1)  electronic  resonance  at  the  D  =  2.87  GHz  zero-field  splitting.  If  a  static  field  is  applied 
with  a  projection  Bn  along  the  NV  axis,  then  to  zeroth  order  the  |  ±  1)  states  are  split  by 
27jvy-B||,  where  ■jnv  is  the  NV  gyromagnetic  ratio  (2.8  MHz/G).  CAW-ODMR  then  records 
resonances  at  E  =  £)±7eRj|,  Furthermore,  if  the  applied  MW  power  is  low  (i.e.  in  the  absence 
of  power  broadening),  the  hyperfine  interaction  with  the  14N  nucleus  of  the  NV  center  can  be 
observed.  The  lineshape  of  the  resonance  in  the  measured  fluorescence  data  depends  on  the 
intrinsic  environment  of  the  diamond  and  the  applied  excitation.  A  single  spin  coupled  to  a 
13C  nuclear  spin  bath  (valid  for  our  low-nitrogen  diamonds)  gives  a  Gaussian  lineshape,  but 
power  broadening  from  the  MW  and  laser  yield  a  Lorentzian  lineshape  The  intrinsic 
linewidth  A v  of  the  Lorentzian  curve  is  related  to  the  inhomogeneous  spin-dephasing  time 
T*  by  A^  =  1/(t rT*). 

For  an  NV  ensemble  in  diamond,  all  four  NV  orientations  are  equally  probable,  giving 
rise  to  as  many  as  eight  electronic  resonances  in  a  CAW-ODMR  spectrum.  However,  if 
the  field  is  perfectly  aligned  to  an  NV  class,  the  projection  of  the  field  along  the  other 


three  classes  will  be  the  same,  leading  to  a  three-fold  degeneracy.  Figure  C-3  (bottom) 


shows  a  measured  CAW-ODMR  spectrum  where  the  field  is  aligned  to  a  single  NV  class. 
In  this  case,  the  aligned  NV  resonances  are  the  outermost  left  (|  —  1)  B  |0))  and  right 
(|1)  B  |0))  transitions  due  to  the  maximal  field  projection  along  the  NV  axis,  while  the 
threefold  degenerate  resonances  from  the  other  three  classes  are  less  split  due  to  a  smaller 
field  projection. 


C.2.4  Laser  Polarization- Assisted  Field  Alignment 

To  assist  in  alignment,  we  take  advantage  of  the  (100)  face  geometry  of  the  CVD  diamond 
and  the  linear  polarization  of  the  exciting  532  nm  laser  beam.  The  NV-  ground  state  is 
preferentially  excited  when  the  laser  polarization  is  perpendicular  to  its  symmetry  axis  [13J. 
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Figure  C-3:  Continuous- Wave  Optically  Detected  Magnetic  Resonance,  (top)  Pulse  sequence 
to  obtain  CAW-ODMR  data,  (bottom)  Measured  CAW-ODMR  spectrum  of  the  diamond 
Noddack  (5  x  1016  e'/cm2,  1250°C  anneal).  The  signal  mi /mo  shows  the  deviation  from 
maximum  fluorescence  when  the  MW  drives  an  NV  orientation  at  or  near  resonance.  This 
data  was  taken  when  the  static  field  was  aligned  to  one  NV  class  (outermost  resonance  pair) 
and  triply  degenerate  with  the  other  three  classes  (innermost  pair).  The  applied  MW  power 
is  set  low  to  observe  the  hyperfine  splitting  of  each  resonance  due  to  interaction  with  the 
14N  nucleus. 
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Figure  C-4:  The  effect  of  laser  polarization  on  NV  class  excitation.  At  one  angle  (130°  with 
the  axes  set  by  the  third  HWP  in  the  setup),  the  green-colored  pair  is  maximally  excited  as 
the  electric  field  is  perpendicular  to  the  NV  axis.  Rotating  the  HWP  45°  in  either  direction 
maximally  excites  the  other  NV  pair.  Selectively  exciting  one  pair  over  another  aids  in 
magnetic  field  alignment. 


Since  the  laser  has  normal  incidence  on  the  diamond  face,  the  polarization  is  in  the  same 
plane  as  the  diamond  face.  |Figure  CV4|  shows  that  when  the  electric  field  of  the  laser  beam 
is  in  the  orientation  shown  on  the  left,  it  preferentially  excites  one  pair  of  NV  classes  since 
it  is  perpendicular  the  NV  axis,  while  also  minimally  exciting  the  other  pair  of  axes  since  it 
is  mostly  parallel  to  those  axes.  Rotating  the  third  half-wave  plate  (HWP)  (just  before  the 
dichroic  mirror)  in  the  experimental  setup  changes  the  laser  polarization,  which  preferentially 


excites  the  other  pair.  |Figure  C-5|  shows  the  CAW-ODMR  spectrum  at  four  different  HWP 
angles.  Around  850°C,  one  pair  of  NV  classes  is  preferentially  excited  (shown  by  the  large 
deviation  in  fluorescence),  while  around  130°,  the  other  pair  is  preferentially  excited.  Going 
left  to  right  in  the  figure,  we  denote  the  first  and  third  NV  resonance  as  pair  1  (preferentially 
excited  at  85°),  and  the  second  and  fourth  as  pair  two  (preferentially  excited  at  130°). 


To  align  the  held,  first  the  NV  classes  in  pair  2  are  made  degenerate  with  the  static  held 
by  means  of  the  translation  stage.  The  resonances  of  both  NV  classes  are  recorded  by  htting 
the  resonance  lineshapes  to  Lorentzians  at  four  positions  of  the  magnet  on  one  axis  of  its 
translation  stage.  The  intersection  of  these  two  lines  yields  the  axis  position  that  renders  the 
two  classes  degenerate.  Once  the  CW-ODMR  signals  from  these  two  classes  coincide  (i.e. 
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Figure  C-5:  CAW-ODMR  spectrum  of  all  four  |0)  -B-  |  —  1)  transitions  for  the  diamond 
Hilbert  (5  x  1015  e  /cm2,  1250°C  anneal)  for  different  half-wave  plate  angles.  The  magnetic 
held  is  mis-aligned  here  such  that  all  four  resonances  are  distinguishable.  Inset:  maximum 
deviation  from  the  reference  fluorescence  of  the  left-most  NV  class  as  a  function  of  the 
half-wave  plate  angle. 


when  the  energy  levels  of  the  two  NV  classes  are  degenerate),  the  HWP  is  rotated  45°  to 
maximally  excite  pair  1.  The  NV  class  in  this  pair  that  we  are  not  aligning  the  magnetic  held 
to  is  now  made  degenerate  with  the  two  NV  classes  in  pair  2.  When  the  resonances  from  the 
three  non-desired  NV  classes  overlap  the  held  is  aligned  and  anisotropic  hyperhne  coupling 
with  the  13C  nuclei  is  eliminated.  It  should  be  noted  that  if  the  diamond  were  isotopically 
enhanced  to  remove  the  13C  content,  the  alignment  of  the  held  would  not  be  necessary  for 
diamonds  with  14N. 


C.3  Resonant  Excitation  and  Field  Projection 


With  the  static  magnetic  held  aligned  to  one  NV  class,  we  now  only  address  this  one  orienta¬ 
tion.  To  implement  the  measurement  protocols  for  T^T^,  and  Tf,  we  now  measure  both  the 
resonance  at  which  the  huorescence  contrast  is  a  maximum  (to  maximize  signal-to-noise)  and 
what  MW  pulse  durations  create  specihc  NV  spin  states.  To  accomplish  the  hrst  task,  we 
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Figure  C-6:  CAW-ODMR  spectra  of  the  (left)  |0)  -B-  |  —  1)  and  (right)  |0)  -H-  |1)  electronic 
transitions  for  the  diamond  Ramon  y  Cajal  (1.6  x  1018  e_/cm2,  850°C  anneal)  The  three 
Lorentzian  peaks  arise  from  the  hyperhne  interaction  with  the  14N  nucleus,  which  can  only 
be  resolved  when  the  NV  ensemble  is  not  power-broadened  from  the  MW  and  optical  drive. 
The  slight  polarization  of  the  peaks  closest  to  the  zero-held  splitting  of  2.87  GHz  arises  from 
partial  polarization  of  the  nuclear  spin. 


take  a  high-resolution  CAW-ODMR  spectrum  around  the  aligned  NY’s  resonance  for  both 


electronic  transitions.  |Figure  C76]  shows  CAW-ODMR  spectra  for  both  of  these  transitions, 
where  each  spectrum  is  fitted  to  a  linear  combination  of  three  Lorentzians: 


g(v)  —  A0  +  (— 


Ai 

/  -  dsy  +  (r/2)2 


A2 

(y  -  /)2  +  (r/2)2 


7I3 

(zz  -  /  +  df)2  +  (r/2)2  ■ 


(C.l) 


Here  g{y)  is  the  normalized  CAW-ODMR  signal,  the  Ai  are  amplitude  ht  parameters,  v  =  / 
is  the  resonance  frequency  of  the  central  resonance,  and  df  is  the  spacing  between  adjacent 
Lorentzians,  usually  2  MHz  for  14N  hyperhne-coupling  (for  15N  only  two  Lorentzians  ht  the 
CAW-ODMR  dip  due  to  the  spin-1/2  nature  of  the  nucleus  [82j).  Both  the  lower  resonance 
zz_  and  higher  resonance  v+  are  extracted,  and  from  the  Zeeman  formula,  the  magnitude  of 
the  static  held  projection  Bn  is  calculated  to  zeroth  order  using 


27  NV 

where  7 nv  =2.8  MHz/G  is  the  NV  gyromagnetic  ratio. 


(C.2) 
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C.4  Rabi  Oscillations 


The  pulse  sequences  used  to  measure  the  coherence  time  scales  require  different  MW  control 
pulses  that  need  to  be  calibrated  for  duration.  For  measuring  Ti,  for  example,  a  MW  n  pulse 
is  needed  to  rotate  the  NV  spin  from  |0)  to  |  —  1).  To  measure  the  duration  of  this  7r  pulse 
(and  other  control  pulses  for  different  measurement  protocols),  we  fix  the  MW  power  and 
sweep  the  duration  of  an  applied  MW  pulse.  Figure  C-7  (top)|shows  the  pulse  sequence  used 
to  measure  population  oscillations:  after  optical  pumping,  a  resonant  MW  pulse  of  fixed 
power  is  applied  for  a  variable  time  r.  During  this  time,  the  NV  spin  population  precesses 
between  the  |0)  state  (low  contrast)  and  |  —  1)  state  (high  contrast).  A  readout  laser  pulse 
measures  the  fluorescence  m i.  The  sequence  is  then  repeated  with  no  MW  pulse  to  provide 
the  reference  signal  tuq.  The  fluorescence  contrast,  defined  for  this  pulse  sequence  as  m\ /mo, 
is  then  plotted  against  the  MW  pulse  duration.  The  signal  is  eventually  damped  due  to 
inhomogeneous  broadening  from  coupling  to  the  environment  and  from  gradients  in  the  MW 
Figure  C-7  (bottom) |  shows  an  example  Rabi  measurement  for  the  diamond  Bohr 


field 


(5  x  1017  e  /cm2,  1250°C  anneal).  The  shape  of  the  curve  is  modelled  as  a  sinusoidal  function 
with  a  stretched  exponential  envelope: 


g{r)  =  A0  +  A\e  {-T,Td)V  cos  (  —  (r  -  r0)  )  , 


2tt 


(C.3) 


where  g(r)  is  the  fluorescence  contrast,  the  A*  are  amplitude  and  offset  constants,  p  is  an 
exponential  stretch  factor,  T  is  the  period  of  the  oscillation,  is  a  decay  time  constant, 
and  To  is  a  phase  shift.  Theoretically,  no  phase  shift  is  required  for  Rabi  nutations;  in  this 
experiment,  however,  minor  technical  issues  with  the  electronics  resulted  in  an  effective  phase 
shift  that  was  constant  for  all  Rabi  nutations  measured  for  all  diamonds. 

To  measure  coherence  times  in  this  thesis  we  required  three  pulse  times:  a  7t/2  pulse  that 
creates  an  equal  superposition  of  the  |0)  and  |  —  1)  spin  states,  a  n  pulse  that  inverts  the 
spin  state  from  |0)  to  |  —  1),  and  a  3ir/2  pulse  that  creates  another  equal  superposition  state 
but  with  a  different  relative  phase.  These  these  times  correspond  to  when  the  cosine  has  a 
phase  of  7r/2,  n,  and  3tt/2  and  we  extract  the  times  from  the  above  fit.  For  the  fluorescence 


signal  in  Figure  C-7(b),  the  times  are  found  to  be  T^/ 2  =12  ns,  Tn  =37  ns,  and  T37r/2  =  61 
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ns. 


C.5  Polarization  Time  Optimization 


We  now  turn  to  optimize  the  duration  of  optical  pumping  by  the  green  laser.  If  the  NV 
electronic  spin  is  in  the  |0)  state  when  excited  by  green  light,  it  will  fluoresce.  If  the  spin 
is  in  |  d=  1),  it  will  initially  fluoresce  less  brightly,  and  any  superposition  of  these  two  states 
will  have  an  intermediate  fluorescence  m  Eventually  the  NV  will  be  pumped  and  these 
signals  converge  into  a  steady  state  value.  The  duration  of  each  applied  green  laser  pulse  is 


chosen  to  always  repump  the  NV  ensemble  back  to  the  |0)  state.  Figure  C-8  (top)  shows  the 
pulse  sequence  used  to  determine  the  minimum  re-pump  time.  After  optically  pumping  the 
NV  with  a  longer  pulse  duration  that  polarizes  the  spin,  a  resonant  MW  ir  pulse  is  applied 
followed  by  the  re-pump  laser  pulse.  After  a  variable  time  r  from  the  start  of  the  re-pump 
pulse  a  fixed-length  gate  pulse  measures  the  fluorescence  m\.  The  sequence  is  repeated  with 
no  7T  pulse  to  measure  fluorescence  mo-  Figure  C-8  (bottom)|  shows  the  two  fluorescence 
signals  as  a  function  of  the  gate  pulse  start  time  r.  The  time  where  these  two  signals  achieve 
the  same  steady  state  value  is  when  the  NV  electronic  spin  is  maximally  polarized  in  the  |0) 
state;  this  time  serves  as  the  minimum  re-pump  laser  pulse  duration. 
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Figure  C-7:  Rabi  nutation  experiment,  (top)  Pulse  sequence  used  to  measure  population 
oscillations  under  a  resonantly  applied  MW  field,  (bottom)  Fluorescence  contrast  mi/m0  as 
a  function  of  MW  pulse  duration  for  the  diamond  Bohr  (5  x  1017  e“/cm2,  1250°C  anneal). 
The  coherent  population  oscillations  yield  a  sinusoidal  fluorescence  signal  with  a  stretched 
exponential  decay  envelope. 
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Figure  C-8:  Optical  pumping  pulse  length  determination  (top)  Pulse  sequence  used  to  mea¬ 
sure  the  re-pump  time  necessary  to  optically  polarize  the  NV  into  the  |0)  state,  (bottom) 
Example  signal  from  a  pump  measurement  for  the  diamond  Lorentz  (5  x  1018  e_/cm2,  1250°C 
anneal).  After  optical  pumping,  a  n  pulse  flips  the  spin  population  and  a  repump  laser  pulse 
is  applied.  After  a  time  r,  the  fixed-length  gate  pulse  is  turned  on  to  measure  the  fluorescence 
m\.  The  sequence  is  repeated  without  the  7 r  pulse  to  give  the  signal  m o.  After  sweeping  the 
start  time  r,  the  time  after  which  both  signals  achieve  their  steady  state  values  is  chosen  as 
the  minimum  duration  of  the  repump  pulse 
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